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THE STATISTICS OF THE GRAVITATIONAL FIELD ARISING 
FROM A RANDOM DISTRIBUTION OF STARS* 


Il. THE SPEED OF FLUCTUATIONS; DYNAMICAL FRICTION; 
SPATIAL CORRELATIONS 


S. CHANDRASEKHAR AND J. VON NEUMANN 


ABSTRACT 


This paper is devoted principally to a statistical analysis of the speed of fluctuations in the force per 
unit mass, F, acting on a star which is moving with a velocity v with respect to the centroid of the 
near-by stars. The solution to the problem depends on the evaluation of the first and the second mo- 
ments of the rate of change of F for a given value of F. 

The statistical problem has been solved on the assumptions of a uniform Poisson distribution of the 
stars and a spherical distribution of the velocities with respect to the chosen local standard of rest. No 
other restrictions have been introduced; in particular, proper allowance has been made for a distribution 
over the different masses M. 


It is found that 
< veF 
= 


where G denotes the constant of gravitation, M the average mass of the stars, the number of stars per 
unit volume, and B a certain function of |F|/Qy (where Qy is a certain “normal” field strength). It is 
indicated how in consequence of this lack of randomness in the rate of change of F for given F and va 
star may experience dynamical friction (i.e., a systematic tendency to be decelerated in the direction of 
its motion by an amount proportional to v ). 

The various second moments of F have also been calculated and lead to an estimate of the mean life 


of the state F. 
The closely related problem of the correlations in F acting at two very close points in a system con- 
taining a random distribution of stars is also considered. 


1. /ntroduction.—In an earlier paper' we analyzed certain statistical features of the 
fluctuating gravitational force acting at some fixed point in a system containing a ran- 
dom distribution of stars in motion. In this paper we propose to extend this discussion 
to the case of the fluctuations in the force acting on a star which is moving with a definite 
velocity v in an appropriately chosen local standard of rest. The discussion of this case 
requires an essential generalization of the analysis contained in I, since for our present 


* Some of the results contained in this paper were included in and formed part of an essay written by 
S. Chandrasekhar and entitled ‘“New Methods in Stellar Dynamics,”’ for which an A. Cressy Morrison 
Prize was awarded by the New York Academy of Sciences. 


1Ap. J., 95, 489, 1941. This paper will be referred to as I. 
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problem what is relevant is the distribution of the velocities of the other stars relative to 
the one under consideration, and we cannot assume that this has a random character,’ 
for the distribution of the relative velocities V will show a marked asymmetry in that 
there will be a preponderance of negative velocities. More exactly, if the distribution of 
the velocities relative to the chosen standard of rest is characterized by randomness, then 


V=-v. (1) 


This asymmetry in the distribution of the relative velocities has important physical 
consequences. Thus, as we shall show (§ 11), it is as a direct result of this asymmetry 
that a star experiences dynamical friction, or, expressed differently, it suffers from a 
systematic tendency to be decelerated in the direction of its motion by an amount pro- 
portional to | v|. 

A second problem we shall consider in this paper concerns the correlation in the forces 
acting at two very close points; this problem is closely related to the one formulated in 
the preceding paragraph. 

2. The general formula for W (F, f).—Consider a star moving with a velocity v. The 
force F acting on the star per unit mass is given by 


where M ; denotes the mass of a typical field star and 1; its instantaneous position vector 
relative to the star under consideration. Accordingly, the rate of change of F is given by 


where V; denotes the velocity of a field star relative to the one under consideration. 
It is clear that the speed of the fluctuations can be specified in terms of the distribu- 


tion function 
W (F,f), (4) 


which gives the simultaneous probability of a given force F acting and an associated 
rate of change of F of amount f. The general expression for this probability can be readily 
written down following Markoff’s method outlined in I, § 2. We have (cf. I, eqs. [18] 
and [19]) 


1 
7 5 
=a ff (p, dpde, (5) 
|p| =0 |o|=0 
where 
14x R'n/3 

limit} 3 
A = drdVdM (6) 

|Vi<e 


In equations (5) and (6) P and @ are two auxiliary vectors, m denotes the average num- 
ber of stars per unit volume, 


an 
tdVdM =7(V;M)dVdM (8) 


2 We use the word “random” in the sense defined in S. Chandrasekhar, Principles of Stellar Dynamics, 
p. 8, University of Chicago Press, 1942. 


| 

| 

| 

| 

| 

(2) | 

| | 

| 
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gives the probability that a star with a relative velocity in the range (V, V + dV) and 
with a mass between M and M + dM will be found. It should be further noted that in 
writing down equations (5) and (6) we have supposed that the fluctuations in the stellar 
distribution which occur are subject only to the restriction of a constant average density. 


Since 


J +drdVdM =1, (9) 


0<M<om Rk |Vi<o 


we can re-write equation (6) as 


limit 


4x R4n/3 (10) 
XK {1 — | 

We replace formula (10) by 

fF 


limit 


R'n/3 
= 1 — (11) 


The integral which occurs in the foregoing equation is conditionally convergent. It 
should be noted that (among others) the integral | . dr is a triple integral, since 


r is a vector. The origin of the expression should make it plausible that one must 
integrate over the two polar angles # and w of r first and over |r| last. (Indeed, 


this ... dr or, rather, .. |r|? sin ddwddd|r| originates from the 
—a 0 0 0 


R 
or, rather, i I sin &dwddd|r|, of equation (10) with 


R—o.) A justification of this plausible procedure will be given elsewhere, by means 
of complex integration. This point is of importance, because an improper order of 
integrations may lead to incorrect results. 

Now equation (11) can also be written as 


A (Py = (12) 


where 
C(p, 


In the foregoing expression for C(P, 7) we shall introduce 9 as the variable of integration 
instead of r. Since @ and r have the same polar co-ordinates and since || and |r| de- 
termine each other, our earlier remarks concerning the r-integration apply equally to the 
$-integration. We have (cf. I, eqs. [22]-[24]) 


dr = . (14) 


fi 1 — . (13) 


| 

| 
| 
| 
| 
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We can re-write equation (13) as 
co +0 +o 
C(p, = f Ji 1 — (5) 
In equation (15) o- will have to be expressed in terms of @. Thus, 
= ||9/2(0-V) } . (16) 
If we now put 
o = (GM)'”o,, (17) 
o-W can be expressed more conveniently as 


Returning to equation (15), we write it in the form 


+0 
C(P, ©) = dMdVrM??D(p, ©) , (19) 
2 ’ 
where 
+00 
D = f {1 — | (20) 


An alternative form for D(p, o) is 
+00 +00 
D = f (1 | + (1 — (21) 


The first of the two integrals which occur in the foregoing equation is equivalent to 
one we have already evaluated in I (eqs. [55|-[58]). We thus have 


+0 


Equations (5), (12), (18), (19), and (22) formally solve the problem of the determina- 
tion of W(F, f). However, an explicit evaluation of W(F, f) would require a complete 
knowledge of the function A(P, o). But if we are interested only in the moments of f 
for a given F, we need only the behavior of A(p, ©) for |o| — 0 or, according to equa- 
tions (12) and (19), in the behavior of D(P, o) for |o|—> 0. We can therefore expand 


(23) 


which occurs under the integral sign in equation (22) in a power series in o. Retaining 
only the first two terms in this expansion, we obtain 


D (Ps ©) = 485 (24) 5? |p Dy ©) + 7) (24) 


where 


+00 
Di (Ps ©) =i || , (25) 


| 
| 
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and 


= 5 fem (26) 


3. The evaluation of D,(P; «).—Substituting for o- from equation (18) in equation 
+ 


(25), we have 
(@-0;) 
D,(P) ©) =ifere} (o,-V) —3 


To evaluate this integral we shall first choose a Cartesian system of co-ordinates with 
the z-axis in the direction of P. Let the vectors o; and V in this system of co-ordinates be 


= So, $3); V= (Vi, Vo, V3). (28) 


(27) 


Further, let 
ly = (/, m, n) = (sin J cos a, sin J sin w, cos 3) (29) 


be a unit vector in the direction of @. Equation (27) now becomes 


CO nr oF 


D, (Ps ©) =if f fd sin 
0 0 0 


x [ (siVit SoVot53V3) —3 (Ps, Vitm’ se Votn?s3V3t+lm[ Ve] 


> (30) 


The integral over w is readily performed, and we obtain 


+1 
(Ps = wif [ $3V3— (s,Vit So V2) ] (31) 


+[$(si:Vit $2V2) —353V3] 2} ||, 


where we have written ¢ = cos 3. Without altering its value we can clearly replace ¢ by 
—t in the foregoing expression. But this replacement changes 


under the integral sign; taking the arithmetic mean of the two resulting integrands, we 
obtain 


co 
fcos( |p| {[ 3 (si Vit | 


(33) 
+[3(s1Vit 
or, writing x = |p| ||, we have 
D,(P) F) = rif fcos (x!) { [s3V3— (51 Vit | (34) 


Vit —353V3] 2} . 


Both the x- and the ¢-integrations can now be carried out. According to the remarks 
made after equations (11) and (13), we must carry out the integration over ¢ first and x 


| 

| ei into (32) 

| 
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later. (We may note here that actually the reverse order of carrying out the integration 
over x first and ¢ later leads to the value 0; this paradoxical result is, of course, caused by 
the — convergence of the multiple integral.) The ¢-integration gives (cf. I, 
eqs. [137]) 


sin x 
x 


(P18) f } (s3Vs— Vit 


0 


i | 
+ GL siVit 351 — [sin x — cos (35) 
= (si Vit V3) f (x?sin x—3sin x+3-x cos 

0 


Now 


dx 
sin x —3 sin x cos x) 


0 


d (1 
2 sin x —3 sin x +3x cos dx 
fo sin x — 3 sin x + 3x cos Wa (=) dx 
(36) 
X Cos x 


1 fd (sinx 1 
x } 
Hence, combining equations (35) and (36), we have 
Dy (Ps = (5, Vi + 52V2—253V3), 


- where, it will be recalled, 5, s2, s; and Vi, V2, V3 are the components of 9 and V in a 
system of co-ordinates in which the z-axis is in the direction of P. If © = (0, o2, o3) in 
the same system, then, according to equation (17), we can re-write equation (37) as 


D,(P, ©) = (GM) (6, o2V2— 203V3). (38) 
4. The evaluation of Do(P, ©).—According to equations (18) and (26), 


+a 


©) =5 (1-V) —3(1e-V) (39) 


where, as in equation (29), 1g is a unit vector in the direction of @. Using polar co- 
ordinates, we can express D2(P, as 


0) =3f f 


(o,-V) —3(lg-V) |? | 
fe [ 1 (40) 


x || ||, ) 
or, introducing a new variable z, defined as 


z= (41) 


| 
| 
| 
| | 

| | 
| 
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we have 
+1 


+1 
0) =3 |p| (1-V) —3 (1g-V) (19-01) (42) 
0 


After performing the integration over w, we obtain 


Da | (o,-V) — 3 (1g-V) |? stds, (43) 


where we have used the bar over [(0,;- V) — 3(1@- V)(1@-9%;)|* to indicate that the aver- 
aging over w has been carried out. 
In order now to evaluate the integral in equation (43) we first note that, since 


[o,-V — 3(le: V)(1@+9;)|? is an even function in ¢ (cf. eqs. [47] and [49] below), 
D.(P,o) has the alternative form 


The integration over z and ¢ are now best performed by regarding them as complex 
variables and integrating along appropriately chosen contours. Thus, writing equation 
(44) as 


“0 


and choosing for s and / paths of integrations as in the appendix to our first paper, we 
obtain 
~1 
. 


In equation (46) the integration over ¢ has to be carried out along a curve from —1 to 
+1 in the complex ¢-plane, which lies entirely in the domain $¢ > 0 and jé| < 1. In 
order to carry out this integration we must first expand [(o,-V) — 3(1@-V)(1@-01)? 
and average over w. We find that 


[ —3(1g-V) (19-03) |? 
=[ (5: Vit 53V3) — 3. P51 + 

+ + mn 53V2) +nl(53Vit 51 V3) 

= (5:V,+ 52Ve+ $3Vs)?—6 (51 Vit (47) 

+n? V3) 


+ 9 siVit V3+n's5 V3) + m2 (s3 Vit $3 Vit 45152Vi V2) 


| | 

| | 
| 

| 
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where 
1=sin cos w ; m =sin J sin w ; n=cos 0(=?#). (48) 
Since 
(1-9); 
m'=3(1-#)?; 


the integral we have to evaluate is 


I = $3V3) — 6 t 53V3) 
X 03 (si Vit (1 + 53 Vat 53 V2) 
(1-2) 20-324 955 428 (ssVit siVit4 Ve) > (50) 
X 2-92 +3 ssVit SiVE+4 
+4535,V3V,) dt. ) 


We readily verify that the various complex integrals which occur in the foregoing ex- 
pression have the following values: 


+1 +1 \ 
frvat=-204); | 
| 
+1 | 


+1 
1) 


1 


| 
+1 +1 | 
| 


1 


Substituting these values in equation (50) and after some minor rearranging of the 
terms, we obtain 


T=—$(1+i) si(S5Vit4V3—2V3) + s3(5Vi+4Vi— 2V3) 


> (52) 
Finally, combining equations (46) and (52) and remembering that 
Re-*4(14+i) = V2 (53) 


and returning to our original variable © = (01, 02, 03) according to equation (17), we 
find that 


Dz (Ps ©) = (24) 3/2 (GM) of (SVit4V2—2V3) 


| 
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5. The expression for A(P, ©) for |o| —> 0.—Combining equations (24), (38), and (54), 
we have 


D (Py 7) = 38; (2m) |p| 9/2 + (GM) 1 + 02V2— 203V3) ) 
(2m) 5/2 (GM) —!|p|-82[ 203) Vit (401+ 503 — 203) V3 | 


+ (403 — 20} — 203) Vi — 80203V2V3 — 80301V 3V 1+ 20;02V iV 2] 
+0(|o|§). 
Substituting this expression for D(P, o) in equation (19), we obtain 
C (Py 7) = (2a) | p | 3/24 ) 
+ (2x) (503+ 403— 203) + (407+ 208) 
x (403 — 203 — 203) — | 
+ +0(|o|%), 


where we have used bars to indicate that the corresponding quantities have been aver- 
aged with the weight function 7(V; M). 
In all the foregoing equations V = (Vi, V2, V3) represents, of course, the velocity of 
a typical field star relative to the one under consideration. If, now, u and v denote, 
respectively, the velocities of the field star and the star under consideration in the chosen 
standard of rest, then 
V=u-v. (57) 


We shall now introduce the assumption that the distribution of the velocities u among 
the stars is spherical,* i.e., the distribution function ¥(u) has the form 


v(u) (58) 


where W is an arbitrary function of the argument specified and the parameter j (of the 
dimensions of [velocity|~') is allowed to be a function of the mass of the star. This as- 
sumption concerning the distribution of the “peculiar” velocities u implies that our 
probability function 7(V; M) must be expressible as 


7(V;M) (59) 


where x(M) governs the distribution over the different masses. For a function 7 of this 
form we clearly have 


MV.=—My; MPV, = (u=1, 2, 3) 
(pn ,v=1, 2,3, 
Substituting these values in equation (56) and after some minor reductions, we find that 
C (Py = (2m) 37 | p | 3/2 — 2 IGM (010, + o202— 20305) ) 
+3 (29) |u|? |p (ai + 02) (24) | 
X [of +403— 203) +03 (40; +5 03 — 203) +03 (403 — 20; — 203) 


8 69030903 — 803010321 + 261020102] +0 ( |o | 3) 


3 It would be entirely feasible at this stage of our work to introduce a more general distribution of the 
velocities (e.g., an ellipsoidal distribution); but we shall not consider these refinements in this paper. 


(61) 


| | 
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where we may recall that (01, o2, 03) and (v%, %, v3) are the components of o and v ina 
system of co-ordinates in which the z-axis is in the direction of P. We shall now specialize 
the co-ordinate system still further by arranging that the vector v lies in the xz-plane 
(see Fig. 1). With this choice of the co-ordinate system 


v,= |v|siny; 72=0; v3= |v|cos7, (62) 


yY=<X (Ps v). (63) 


where 


The expression for C(P, ©) now simplifies to 


C (Py = (24) 32 |p |3/2— | v| (o,siny—203cosy) 


+1(27) 3/2G1/24f 1/72 lu |? | p | —3/2 a3) + (27) 3/2G1/2f172 | v | 2 | p | —3/2 


{ 
(64) 
[of (5 sin? y — 2 cos? y) +03 (4sin?y — 2cos?y) +03 (4 cos? y — 2 sin 
— 80,03 sin y cos y] +O(|o|?). J 
Substituting this expression for C(P,@) in equation (12), which defines A(P,o), we obtain 


where we have written 


1/2 (66) 4 
g=2nGM \u\n; bag 


and 

=a;(5 sin? y — 2 cos y) +03 (4 sin? y — 2 cos? y) (67) 
+ 03 (4 cos? y — 2 sin? y) — 8030, sin y cos y. 

An alternative form for A(P, ©), which we shall find useful, may be noted here: 


+kR(o) J 


Now, according to equation (5), A(P, o) is the six-dimensional Fourier transform of 
the distribution function W(F, f). Consequently, for the purposes of this equation the 
vectors P and o should be referred to a fixed system of co-ordinates. But the expressions 
for P(o), Q(@), and R(o) (eq. [67]) involve the components of o referred to a variable 


‘Tt will be noted that our present definitions of a and 6 agree with those of our earlier paper (I, eq. 


| 
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system of co-ordinates, depending on the direction of P. We shall now give the linear 
transformation required to pass from this variable xyz-system to a fixed énf-system (see 
Fig. 1). This fixed &nf-system is so chosen that the ¢-axis is in the direction of F and 


Fic. 1 


the é{-plane contains the vector v. The accompanying table gives the direction cosines of 
the axes belonging to one system with respect to the axes belonging to the other. 


On OF 


sina —Ilcosy | m COS ¥ cos a — n cos y 
sin y sin y sin y (69) 
m cos a nsina—Icosa | m sin a 
| sin y | siny sin y 
As = sin 9 cosw | =sindsinw = m = = 


Here a is the angle between F and v: 


a=<X(F,v), (70) 
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and 
cos y=cos X(P) v) =ncosat/sina. (71) 


Thus the required linear transformation is 
(¢=1,2,3) (72) 


where g¢, o,,and a; are the components of 9 referred to the fixed system of co-ordinates. 
6. The evaluation of the first moment of f—To determine the average value of f (for a 

given F) in any specified direction, it would clearly be sufficient to evaluate the first 

moments of the components of f (namely, /:, f,, and f¢) along the three principal direc- 

tions of the énf-system defined in § 5. We shall consider first the moment of /:. 
According to equation (5), 


1 +0 +0 


+00 
[WED ff fee fedPdods. (73) 


But 
$00 


1 | 
= 15" (04) 6 (oy) (74) 


where the 6’s denote Dirac’s 6-functions and 6’ the first derivative of the 6-function. 
Remembering that 


+0 
J f(x) (75) 
we can reduce equation (73) to 
+00 +00 
7 = F o | 7 
But, according to equation (68), 
OP 
since P is linear ino. Thus, 
+00 +00 
Jw (F,f) fedf (78) 


or, choosing polar co-ordinates (see Fig. 1), 


feaf anf fe jay) dd |P|dddw. (79) 


Performing the integration over w, we obtain 


+00 +1 aD 
[Wes = 85 f (22) plea |p at, (80) 


| 
| 
| 


| 
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where we have used a bar over 0P/do; to denote that the averaging over w has been 
carried out; further, in equation (80) we have changed from the variable 3 to ¢ = cos #. 


Putting (cf. I, eq. [134]) 
[P| |F| (81) 


in equation (80) and remembering that 0P/da¢ is even in / (cf. eq. [87] below), we obtain 
+00 | aD. 
7 (+ 2 
Jw (PA) fedf = x?cosxtdxdt, 


The first moment of f: is now obtained by dividing the foregoing equation by W(F) 
(1, eq. [117]). We thus have 


fe= [fe ~ (2 /8)*/ x? cos xtdxdl. (83) 


We have similar expressions for f, and fy. 
We shall now evaluate dP/dez, etc., According to equation (67), 


P= o,sin y— 203 cosy, (84) 
or, transforming to the énf-system according to equation (72), we have 
P= (Ayoe+ +r sin y — 2 (Agoz + +307) COS (85) 


Using the table of direction cosines (69) and after some rearranging of the terms, we 
obtain 


P= (1 — 3n?) cos a — 3ln cos sin a) (86) 
86 
—3o,m(lsina+n cosa). 


From this equation we readily find that 


oP = —3(mlsina+mncosa) =0, (89) 
= (1—3n*) cos a—3lnsina=+(1— 3) cosa. 


Thus, according to equations (83) and (87), we have 
(By 4 (8) sin @ fr =0; (8) (88) 


where we have written 


(1 — 322) cos xtdxdt. (89) 


{ 
| 
{ 
| 
q 
| 
| 
| 
| 
| 
( 
{ 
i 
| 
q 
| q 
| 
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We shall now evaluate $(8). Performing the integration over ¢, we obtain 


(B) fe x (2x cos sin x 2sin x) de (90) 


or, after some rearranging of the terms, 


J (B) =6 | (sin x — x cos x) 2 sinxdx. (91) 
0 < 0 


The second of the two integrals which occur in equation (91) is seen to be related very 
simply to the Holtsmark function (8), defined in I, equation (116). Further, denoting 


by A(8) the integral 


dx 
: =~ f (sin x — ) 
K (B) cS (sin x — x cos x) = (92) 
we have 
J(B) —7BH(B). (93) 
Now 
= (sin x — x cos x) | 
: | (94) 
| 
cos x) dx 
) (sin x — x cos v) dx, | 
or, after integrating by parts, we have 
ale rsin = (95) 
Hence 
fi 
K (8) = | A(8)dp. (96) 
0 
Now, combining equations (88), (93), and (96), we have 
fe = —gB(B)sin a ; fr=0; fr = 2gB(B) cosa, (97) 
where 
a8 
(98) 


The asymptotic properties of B() are easily derived from those of H(8) (1, eqs. [118| 
and [119]). We find that 
B(B) = C3") B? +0 (B—0) (99) 


and 
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Consider now an arbitrary direction (J, m, n). Then 


(101) 


or, according to equation (97), 
tiv min = —3 (8) | v| (Lsin a— 2n cos a), (102) 


where we have substituted for g from equation (66). Remembering that the direction 
cosines of v are (sin a, 0, cos a), we have 
(lsina — 2n cosa) =|v| (Jsina+ncos a) —3n|v|cosa 


103) 


where 1),,, Stands for a unit vector in the direction (J, m, n). Hence we can re-write 
equation (102) as 


or, more simply, as 


We shall return to this equation in § 11. 
7. The evaluation of the second moment of f.—There are in general six independent mo- 
ments of the second order that we have to consider, namely, 


In terms of these six moments, the second moment of the resolved component of f along 


any arbitrary direction can be specified. For, if fi,m,, denotes the component of f in the 
direction (/, m, n), then 


Ufet mfatnfr)? (107) 
=PRitm fit 


In our particular problem we should expect on symmetry grounds (and this is verified to 
to be the case) that two of the six moments (eq. [106]), namely, fe fy and Safes vanish 
identically (cf. eq. [97], according to which f, = 0). Accordingly, equation (107) reduces 
in our case to 


=P fit 2nlfrfe. (107) ’ 


m,n 


We shall first consider the moments fi, ra and i We have (cf. I, eq. [123]) 


fw (F, fe A | dp. (108) 


=0 


j= (21) (105) | 
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where we have used 7 to denote either £, n, or £. Now, according to equation (68), 


0° 
| (109) 
Hence, 
+00 
fw (F, f) frdf = 
(110) 
or ant 
or, choosing polar co-ordinates in the (£, n, ¢)-system, 
+120 
(111) 
x 2 OF + |p|. 


Performing the integration over the azimuthal angle w and introducing the variable 
x= =0%6, (112) 


equation (111) becomes 


+ 4abp3/ +e (3 k cos xtdtdx , | 


+00 
JWOD 5 
= (113) 


where the bars over 0?P?/d0?, etc., indicate that the averaging over w has been carried 
out. We shall presently show that 0?P?/d0?, etc., have the forms 


\ 


#0 | 
( = (r= &, Ns (114) 
where Po,;,...., Rs, are constants with respect to the variables of integration in equa- 


tion (113). Substituting these equivalents for 0°P?/do?, etc., in equation (113) and di- 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
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viding throughout by W(F) according to I, equation (117), we obtain for the moment of 
f? the expression 


1/72 
vie (115) 


9 1 
3/2 2 | 
(B) {fe st? x? cos (xt) dxdt, | 
where, for the sake of brevity, we have written 


In equation (115) the integration over ¢ can now be effected. Using the elementary 
formulae 


1 P 
sin x 
feos xtdt = 
A x 
ms 2 
fe cos xtdt = ——,(sinx—x cos x), (117) | 
A x x 
sinx 4 | 
cos xtdt= (—x* cos x + 6x cos x + 3x? sin x —6 sin x) , 
x / 


“0 
we readily obtain from equation (115) 


ab (Ly, Lu )G(B) — (B) J (8) | 


H (Bp) | 
(118) | 
| 
26H (B) (Po, BH (B) — 2P2, -K (B) — 4 Px, (B) | 
where G, H, I, J, K, and L are functions of 6 defined by 
) 
=— sin xd x 
G (8) sin xdx , 
= sin xdx , 
mB 
») @ 
=— (sinx—x cos x)dx, 
0 
oo (119)5 
J(B) = fe ~ (— x3cos x + 6x cosx + 3x*sin x — 6sinx) dx, 
“0 
K = fe (Sinx — x cos x) dx, 
0 


co 
(—xcosx+6xcosx+3x?sin x —6sin x) dx. J 
0 


5 We may draw attention to the fact that we have not changed here our definition of H(8); it repre- 
sents the Holtsmark function as hitherto; similarly, our present definitions of G(8) and K(8) agree with 
those we have given earlier (I, eq. [159] and eq. [92] of this paper). However, the functions (8), J(8), 
and L(8) are introduced here for the first time. 


(| | 
| 
| 
| 
| 
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It remains to obtain the expressions for Lo,;,.... (7 = n, The calcula- a 
tions are straightforward but somewhat tedious. But we shall illustrate the method by g 
outlining the derivation of (Ro,¢, ‘Re,z, and 

According to equation (67) g 


(0;[5 sin? y — 2 cos? y] + sin? y — 2 cos® y] | 

oe, 99; > (120) 
+ o3[4 cos? y — 2 sin? y] — 8030; sin y cos y). } 

To carry out the differentiation we must first apply to (1, 2, 03) the linear transforma- 

tion (72). It is, however, not necessary to carry out this transformation explicitly. For, 

since clearly 


2 
do; do; 
we have 
——= 2[dj(5 sin? y — 2 cos? y) + A3(4 sin? y — 2 cos? y) | 
do; 


+3 (4 cos? y — 2 sin? y) — 8A4A3 sin cos 
or, after some minor reductions, 
2[—2+6)3 cos? y + (7A; + 6A3) sin? y — 8A,A3 Sin y cos y}] (123) 
g 
Substituting for Ay, Xe, and A; from the table of direction cosines (69) and using for cos 
its equivalent (71), we obtain, after some rearranging of the terms, 
—.= 2[-24+ sin? a+ 21Pn? cos? | 
do; (124) 
+ 6m? cos? a — 227 sin? a — 22In sin a cos a], | 


where it will be recalled that a is the angle between F and v. Averaging the quantity on 
the right-hand side of equation (142), we finally obtain 
(2+ sin? a) + (15 —4 sin? a) P+ (—214+44% sin? a), (125) 


do, 


from which the values of ‘Ro, ¢, ‘Re,z, and ‘Ky, can at once be read. The evaluation of 
the other quantities proceeds similarly. Collecting all the results, we obtain the following 
table of values: 


Po, sin? a; Pr», gm +9 38 sin? a: t= + 6,8 sin? a. 

Po sin? a; —O9+ sin? a, (126) 

Po 2cos?a; Po, =+18— 27 sin’ a; 


| 
| 
| 
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and 
Ro ge=2 + isin? a; Rey ¢=t+15 sin? a; Ry ¢=—21+ sin? a 
a; Ro +15 sin? a; Ry —21+ 435 sin? a, (128) 


Ro, =10—8 sin? a; =—444+59sin2 a; =t42— 63 sin? a 


We have still to evaluate the cross-moments ils fafe, and cfg. As we have already 


stated, the first two vanish identically, and we have only to consider f; f¢. Analogous to 
equation (110) we now have 


1 +0 
fu (F, f) fefedf =z; fee (129) 


or, changing to polar co-ordinates and integrating over w, we find (cf. eq. [113]) 


(130) 


| + cos xtdtdx , 


where we have made a further change of variables according to equation (112). In 
equation (130) the bars are again to be understood as indicating that the corresponding 
quantities have been averaged over w. An elementary, but somewhat lengthy, calcula- 
tion yields 


(5 = ( —1+15f—18?')sinacosa, | 
(131) 
= (-3+ 377 — 42#) sin a cosa. | 
Oordor 
Further, 
(132) 


Substituting these formulae in equations (130) and dividing throughout by W(F) ac- 
cording to I, equation (117), we obtain 


2 
rH (B) 


(B) fhe cos xtdtdx{ sin a cos a, 


Sf —42t) cos xtdldx 


(133) 


| — 
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or, integrating over ¢ and introducing the auxiliary functions G,...., L, according to 
equation (119), we have 


[-8G (8) — 741 (8) +168 (8) ] 


H(B) | 
(134) | 
+3597 48H (8) — 30K (8) + 72L(8) ] sin a cos a. 
Equations (118) and (134), together with 
(135) 


represent, then, our expressions for the six independent second-order moments. 


8. The moment |f|* and its average fi? The mean life of the state F.—In § 7 we evalu- 
ated all the independent second-order moments that exist and, as we have already 
stated, the second moment of the resolved component of f along any direction can be 
specified i in terms of them. However, the most important quantity is the average value 
of |f|*. This is clearly given by 


t= > F (136) 


4 


or, according to equation (118), by 


Pi +) L(B) | 


r=E, 0, ¢ 

On the other hand, from equations (126)-(128) we find that 

Po, 5 sin? a+ 2 cos? a ; -=0, 

; 2, -=0 ; =9, (138) 
>) Ro, += 14—7 sin? a : Rs, -=—14+421 sin? a; | 


Again, according to equations (138), 


=8—6sinta; 9,,=4; 14 sin? a. (139) 


T Tt 


We thus have 


(8) +7 (sin? aG(B)—(3 sin? a—2)/(B) 
(140) \ 


{ 
| 
| 
| 


(4—3 sin? a) BH (8) +3(3 sin? a — 2) K(8) 


\ 
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The foregoing expression gives the mean square value of the rate of change, f, to be 
expected in the intensity of the force, F, acting on a star when it is further known that 
the direction of F makes an angle a with the direction of motion. But the possible direc- 
tions of F, for a given direction of motion, are distributed uniformly over the unit sphere. 
Hence, the mean square value of f for a given value of F and for all relative orientations 
of the two vectors F and uv is obtained by simply averaging equation (140) over a. We 


thus obtain 


or, substituting for k and g?/2ab from equation (66), we find that 


fl? jl = 448 H (1 5) +35 (142) 


If |v; — 0, we recover the formula of our earlier paper (cf. I, eq. [158]). 
The formula (142) has an immediate application for estimating the mean life of the 
state of fluctuation in which a force F per unit mass acts on a star moving with a speed 
v.. For, arguing as in I, § 9, we may define the mean life by the equation 


(143) 


According to equations (112) and (142), we therefore have 


(144) 
_H (6) 
(B) 
From this equation we derive the relation 


9. The auxiliary functions G, H, 1, J, K, and L.—In §§ 6, 7, and 8 we have seen that 
our expressions for the various moments involve some or all of the functions G(8), H(8), 
1(8), J(8), K(8), and L(8). These functions are all defined by certain definite integrals 
(eq. [119]), the integrands of which contain 8 as a parameter. We shall now show how 
the functions G(8), 7(8), J(8), A(8), and L(8) can all be expressed in terms of the 
Holtsmark function, //(8). We shall also establish certain relations which exist among 
them. 

We have already seen that the functions G(8) and A(8) are related to H(8), according 
to the formulae (eq. [96] and I, eq. [162]), 


G(B) = fe 3/211 (B) dB (146) 
and 
B 
K (B) = a8. (147) 


It remains to establish similar relations for 7, J, and L. 


| 
(141) | 
| 
\ 
) 
\ 
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Consider, first, 7(8). Writing the integral, defining it in the form 
=-=— fe {x/8)*/* (sin x — x cos x) q 
dx 


and integrating by parts, we find 
4 fee) 


(B) 


3m 0 


or, remembering our definition of G(8) (eq. [119]), we have 


co 


(B) =3G(g) — 


B3/2 
But 
fe (sin  — x cos x) . 
Hence, 
dl 


in other words, /({) satisfies the differential equation 
—+37=G(B). 
B dp +3 8) 
The solution of equation (153) appropriate for us is 


1(8) (B) dB. 


e~ (7/6)? y—1/2 sin xdx—— _ [| e~ (sin x — x cos x) 
0 


fe (x/6)*/* (sin x —x cos x) 


(148) 


ix 


(149) 


(150) 


This formula is useful for the purposes of numerically evaluating the function /(£). 
We shall now establish for /(8) a relation similar to equation (154). We have 


J(B) =—-— [Je7'/8*? (—x3 cos x + 6x cos x + 3x? sin x — 6 sin x) 
in 
1 
“_(x-7/2) dx 
dx 
= sin xd x 
iw A ( 1 5) 
aa fe (— x3 cos x + 6x cos + 3x? sin x — 6 sin ds 
/ x3 | 
dJ 
= 2G (f 28 —., 
7(B) 8 dp 
Hence /() satisfies the differential equation 
= Gp) (156) 


dg 


| 

\ 

| 
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Consequently, 
B 
0 


This is the required formula. 
By following a procedure exactly similar to that adopted for treating /(8) and J(6) 
we can show that L(£) satisfies the differential equation 


- = BH (B) (158) 
and that therefore 
> B 
(B) dB. (159) 
0 
Further, according to equations (150) and (155), we have the relations 
I (B) = 3G(B) —B-*”K (B) (160) 
anc 
J (B) = 7G(B) (161) 
Finally, we may note the following asymptotic forms for the various functions: 
| 4 2 
G(B) 3_ G (8) 
4 
H(8) > H (B) 6 
162 
K (8) > 1 
4 
Li) L(g) > 
oT 
B (B) ( 1,0) B(B) 3 B3/2 


10. The correlations in F acting at two very close points.—The formal theory developed 
in the preceding sections has direct applications to a different problem, namely, that of 
the correlations in the force acting at two very close points, for the difference between 
the values of F acting at two points distant 6r from each other is given by 
ér 6r;) 


6 
rac 


AF = c> 


where we have assumed that one of the points is at the origin of our system of co- 
ordinates. Comparing equations (3) and (163), we see that formally the problems of 
specifying the distributions W(F, f) and W(F, AF) differ only to the extent that, while 
in the first case we have to allow for a distribution over the relative velocities V, in our 
present problem 6r is a fixed constant vector. 


§ We discuss later in this section the validity of this formula. 
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Thus, expressing W(F, AF) in the form 


+c +c 


W (F, AF) AF) (py , (164) 


we have (cf. eq. [12]) 


where (cf. eq. [15]) 


0 


In equation (166) W now stands for (cf. eq. [16]) 
From this point on, the analysis proceds exactly as in §§ 3, 4, and 5. It is thus seen that 
equation (64) is now replaced by (cf. eq. [56]) 
C (Py %) = | 4+ 2 | sin y — 203 cos y) 
+ (29) | or |? |p| [67 (5 sin? y — 2 cos? y) +03(4 sin? y > (168) 


— 2 cos? y) + 03 (4 cos? y — 2 sin* y) — 80,03 sin y cos y] +O( , 


where 
Y=<X (P, or) (169) 


and the co-ordinate system has been so chosen that the z-axis is in the direction of P and 
6r lies in the xz-plane (see Fig. 1). Hence, 


A (Py ©) = ipl + (170) 


Ww here a, P(o), and R(o) have the same meanings as in equations (66) and (67), while 0’ 
and g’ now os for 


The evaluation of the first and the second moments of AF proceeds as in §§ 6, 7, and 8. 
Thus, analogous to equation (105), we now have 


AF r, 3xGMnB (5, )(#-3 F). (172) 


Similarly the moment of | AF |? is given by (cf. eq. [140]) 


= 14ab’ aG (8) — (3 sin? a | 


| 
| 
| 
| 
| 
, | 
(173) 
+ ey “4 3 sin? a) BH (B) +3(3 sin? a 2) K(B)], | 
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or, substituting for a, 6’, and g’ from rt (66) and (171), we have 


= | or 


5 [sin? aG (B) 


[ (4—3 sin? a) 


GM | ar| 
— (3 sin?a 2) 1(8) 144% 


x BH (B) +3(3 sin? a—2)K(B)].) 


Averaging the foregoing expression over all a, we find 


M | 6r |? 


Fi, 
Now, according to equation (172) for a fixed |6r|, AF tends to infinity as | F| > ~. 


This is contrary to what we should expect on physical grounds, namely, that AF should 
tend to zeroas | F| — , for, since the highest fields are approximately produced by the 
nearest neighbor,’ it follows that, as |F| — ©, the particular star which effectively 
produces the field must be so close to one of the two points considered that no correla- 


tions in the directions of F acting at the two points can be expected. In other words, AF 
should tend to vanish as | F| — ~. But this same argument shows why our present 
theory of spatial correlations fails as | F| — ©, for, given a | 6r|, however small, we can 
always choose a |F) so large that on the first-neighbor approximation, the “nearest 
neighbor’ will be closer to one of the points than | ér|. Under these circumstances the 
contribution to AF arising from this nearest neighbor will no longer be represented to any 
degree of accuracy by a term of the series (163)— the Taylor expansion of r/|r|* on 
which this series is based will cease to be valid for at least that particular term cor- 
responding to the nearest neighbor producing a |F| — . We therefore conclude that 
our present method gives only the asy mptotic hehawies of the true spatial correlations 
in the sense that, given a | F|, however large, we can choose a | 6r| sufficiently small for 
our formulae to be valid for | | Fl less than the specified limit. 

11. Dynamical friction —We ‘shall now return to equation (105) for a discussion of 
its implications for general dynamical theory. According to this eins 


— 


where B(8) is defined in equation (98). We shall first derive certain formal consequences 


of this equation. 
Multiplying equation (176) scalarly with F, we obtain 


+ or|?. (175) 


dF\ (‘FL 
But 
dF 
Hence 


7S. Chandrasekhar, Ap. J., 94, 511, 1941 (§§ 3 and 4). 


) > 
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On the other hand, if /; denotes the component of F in an arbitrary direction at right 
angles to the direction of v, then, according to equation (176), 


ar; 180 
27GM nB On ( ) 
Combining equations (179) and (180), we have 

1 ‘ee | 

F;\ dt/r,v 2 |F| 
Equation (181) is clearly equivalent to 

| 

[5 = 0. (182) 


In other words, we have proved that 


=(). 183 
dt |F | 3/2 F, 


We shall now examine the physical consequences of equation (176) more closely. In 
words, the meaning of this equation is that the eg of 


(5+) Tae (184) 


along any particular direction gives the average value of the rate of change of the force 
per unit mass acting on a star that is to be expected in the specified direction when the 
star is moving with a velocity v in an appropriately chosen local standard of rest. Stated 
in this manner, the essential difference is at once seen in the stochastic variations of F 


with time in the two cases |v| = O and |v| # 0. In the former case F = 0; but this is 
not generally true when |v, + 0. Or, expressed somewhat differently, when |v) = 0 
the changes in F occur with equal probability in all directions, while this is not the case 
when |v, # 0. The exact nature of this difference is brought out quite clearly when we 
consider 


8: 
( dt (285) 


according to equation (179). Remembering that B(8) 2 0 for 8 > 0, we conclude from 
equation (179) that 


(186) 
dt F,vu 
and 


In other words, if F has a positive component in the direction of v, F increases on the 
average; while, if F has a negative component in the direction of v, F decreases on the 
average. We shall now show that it is this essential asymmetry introduced by the direc- 
tion of v that gives rise to the phenomenon of dynamical friction. 


| 
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The characteristic aspects of the situation governed by equation (179) are best under- 


stood when we contrast it with the case F = 0. Under these circumstances we can 
visualize the motion of the representative point in the velocity space somewhat as fol- 
lows:’ The representative point suffers small random displacements in a manner that 
can be adequately described by the theory of random flights.* More specifically, the star 
may be assumed to suffer a large number of discrete increases in velocity of amounts 
F1T(F), where 7 is the ‘‘mean life’’ of the state F; these increases are further assumed to 
take place in random directions. On these assumptions it readily follows that the mean 
square increase in the velocity, which we may expect the star to suffer in time ¢, is given 
by 


|F[°7¢. (188) 
An alternative way of describing the same situation is to assert that the function P(y, é; 


vo), which gives the probability that the star has a velocity v at time #, given that v = uo 
at ¢ = 0, satisfies the diffusion ae 


Gat av? (339) 
with the “coefficient of diffusion,” D, having the value 
(190) 
The solution of equation (189) appropriate for our purposes is 
1 
t Uo) (4x Dip / (191) 


The formula (188) is seen to follow immediately from the foregoing solution. 
Returning to the discussion of the case governed by equations (176) and (179), we 
see at once that the idealization of the motion of the representative point in the velocity 
space as a true problem in random flights can no longer be valid. Moreover, according 
to equations (183), (186), and (187), during a given state of fluctuation the star is likely to 
suffer a greater amount of acceleration in the direction of — v when (v- F) is negative than 
in the direction of +v when (v- F) is positive. But the a priori probabilities for (v- F) 
to be positive or negative are equal. Hence, when integrated over a large number of fluc- 
tuations, the star must suffer cumulatively a larger absolute amount of acceleration in 
the direction opposite to its direction of motion than in the direction of motion. In 
other words, there is a tendency for the star to be relatively decelerated in the direc- 
tion of its motion; further, this tendency is proportional to | v|. But these are exactly 
what are implied by the existence of dynamical friction.'® 
YERKES OBSERVATORY 
WILLIAMS Bay, WISCONSIN 
AND 
INSTITUTE FOR ADVANCED STUDY 
PRINCETON, NEW JERSEY 


8 Cf. ibid., $§ 2 and 7. * See, e.g., Lord Rayleigh, Phil. Mag., 6th ser., 37, 321, 1919. 
1° See two forthcoming papers by one of us (S.C.) on ‘Dynamical friction” in an early issue of the 
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NUCLEAR EMISSION IN SPIRAL NEBULAE* 
K. SEYFERTT 


ABSTRACT 


Spectrograms of dispersion 37-200 A/mm have been obtained of six extragalactic nebulae with high- 
excitation nuclear emission lines superposed on a normal G-type spectrum. All the stronger emission 
lines from 4 3727 to \ 6731 found in planetaries like NGC 7027 appear in the spectra of the two brightest 
spirals observed, NGC 1068 and NGC 4151. 

Apparent relative intensities of the emission lines in the six spirals were reduced to true relative in- 
tensities. Color temperatures of the continua of each spiral were determined for this purpose. 

The observed relative intensities of the emission lines exhibit large variations from nebula to nebula. 
Profiles of the emission lines show that all the lines are broadened, presumably by Doppler motion, by 
amounts varying up to 8500 km/sec for the total width of the hydrogen lines in NGC 3516 and NGC 7469. 
The hydrogen lines in NGC 4151 have relatively narrow cores with wide wings, 7500 km/sec in total 
breadth. Similar wings are found for the Balmer lines in NGC 7469. The lines of the other ions show no 
evidence of wide wings. Some of the lines exhibit strong asymmetries, usually in the sense that the violet 
side of the line is stronger than the red. 

In NGC 7469 the absorption K line of Ca 1 is shallow and 50 A wide, at least twice as wide as in 
normal spirals. 

Absorption minima are found in six of the stronger emission lines in NGC 1068, in one line in NGC 
4151, and one in NGC 7469. Evidence from measures of wave length and equivalent widths suggests that 
these absorption minima arise from the G-type spectra on which the emissions are superposed. 

The maximum width of the Balmer emission lines seems to increase with the absolute magnitude of 
the nucleus and with the ratio of the light in the nucleus to the total light of the nebula. The emission 
lines in the brightest diffuse nebulae in other extragalactic objects do not appear to have wide emission 
lines similar to those found in the nuclei of emission spirals. 


Many of the spectra of extragalactic nebulae obtained at the Mount Wilson and Lick 
observatories show one or more emission lines in addition to the usual absorption spectra. 
In particular, N. U. Mayall! finds that 50 per cent of his spectra of spirals show the 
[O 11] doublet \ 3727 in emission either in the nuclear region or in the arms. However, 
only a very small proportion of extragalactic nebulae show spectra having many high-ex- 
citation emission lines localized in the nuclei. These emission features are similgr to those 
found in planetary nebulae and are superposed on the characteristic solar-type absorption 
spectra. Twelve nebulae? are now known which probably belong to this unusual class of 
objects. Most of them are intermediate-type spirals with ill-defined amorphous arms, 
their most consistent characteristic being an exceedingly luminous stellar or semistellar 
nucleus which contains a relatively large percentage of the total light of the system. 
Plate I shows a photograph of NGC 4151, a typical example of this type of nebula. 

Probably the earliest spectrographic observation of a member of this unusual class of 
objects was that of NGC 1068 by E. A. Fath,’ in which he found five emission and two 
absorption lines. In 1917, V. M. Slipher* found hydrogen lines and the nebular lines 
N1 and N2 bright in the nucleus of NGC 5236. Shortly afterwards, Slipher’ made, the 
discovery that the emission lines in NGC 1068 were not monochromatic images of the 
. slit but were small ‘‘discs.”’ These findings were confirmed by Campbell and Moore,* 


* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washington, No. 671. 
+ Fellow of the National Research Council. 

1 Lick Obs. Bull., 19, 33, 1939. 

2 NGC 1068, 1275, 2782, 3077, 3227, 3516, 4051, 4151, 4258, 5548, 6814, and 7469. 

3 Lick Obs. Bull., 5, 71, 1908. 

4 Pop. Astr., 25, 36, 1917; Proc. Amer. Phil. Soc., 56, 403, 1917. 

> Lowell Obs. Bull., 3, 59, 1917. 6 Lick Obs. Pub., 13, 88, 1918. 
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who stated that the bright bands in NGC 1068 were “‘fully 30 Angstroms wide.’’ Camp- 
bell and Moore’ also found in NGC 4151 emission lines several angstroms wide. E. P. 
Hubble® refers to planetary-type emission in the nucleus of NGC 4051 (as well as 
NGC 1068 and NGC 4151) in a paper published in 1926. Detailed descriptions of the 
typical emission spirals NGC 1275 and NGC 4151 were published by M. L. Humason® 
and N. U. Mayall,'® respectively. In addition, a private communication from Dr. Mayall 
indicates that NGC 3516 and NGC 7469 show broad bands of hydrogen in emission and 
hence belong in the group under consideration. 


I. THE OBSERVATIONAL MATERIAL 


The present investigation is an intensive study of six of the brightest extragalactic 
nebulae showing emission bands in their nuclei (Table 1). Of these six, special emphasis 
was placed on the three having the brightest nuclei, NGC 1068, 3516, and 4151, because 


TABLE 1* 
EMISSION SPIRALS OBSERVED 


1950 

NGC TYPE otal | SPECT. | Mopuzus 

R.A. Dec. | | 
1068............| 2640.1 | — 0°44 Sb 100 13.0 G3 | 2670 17 
T on 315.6 | +41 18 E: 13.0 | 15.5 G3 | 30.0 4 
3516... | | Sa 12.2 | 13.7 G2: | 28.5 6 
4051... .. 12 0.6 | +44 48 Sb 11.7 140 G2 | 260 4 
4151... 12 80 | +39 41 Sb 112 120 | G2 | 26.0 12 
..... 3-67 | Sa 13.0  14.3:| GO: | 29.8 2 

| | 


* The total apparent photographic magnitudes are from the Shapley-A mes Catalogue of External Galaxies (Harv. Ann., 88, 43, 
1932). The apparent magnitudes (photographic) of the nuclei were estimated from short-exposure plates, taken in series with 
selected areas. The distance moduli are new determinations derived from magnitudes of resolved stars in the arms (NGC 1068), 
radial velocity (NGC 1068, 3516, 7469), or from association with recognized clusters or groups (NGC 1275, 4051, 4151). The 
plates used for determinations of nuclear magnitudes and most of the data for computing the distance moduli were supplied by 
E. P. Hubble. The spectral types were determined by M. L. Humason. 


it was possible to observe them with higher dispersion than could be used on the fainter 
objects. 

The observations were obtained mainly at the Cassegrain focus of the 60-inch re- 
flector using (1) the one-prism spectrograph and cameras of 3.5, 6, 9, 10, and 18 inches in 
focal length, and (2) the two-prism nebular spectrograph ‘‘VIA” with a 3-inch camera. 
In addition, a few plates were obtained with a grating spectrograph and Schmidt camera 
and with the “VIC” Cassegrain two-prism spectrograph and 10-inch camera at the 
100-inch reflector. Table 2 lists the dispersions with the various combinations. 

Spectral calibration strips for determining relative intensities at each wave length 
were placed beside the nebular spectrum on the same photographic plate during the 
course of the nebular exposure. These standardization exposures were never shorter than 
5 minutes and were usually 30 minutes in length, or 13 per cent as long as the average 
nebular exposure. For a few plates the ratio of standard to nebular exposure was as low 
as 1 or 2 per cent. Examination showed that the photometric results obtained from 
these plates were not systematically different from plates with much larger ratios. 


7 Lick Obs. Bull., 13, 122, 1918. 
5 Mt. W. Contr., No. 324; Ap. J., 64, 328, 1926. 
9 Pub. A.S.P., 44, 267, 1932. 10 Pub. A.S.P., 46, 134, 1934. 
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Hence we believe that the relatively short standardization exposures do not introduce 
any serious errors. 

The chief photographic plates used were Cramer Hi-Speed, Eastman 103-O, and 
Eastman 103-F. To increase their sensitivity, the two blue-sensitive emulsions (Hi- 
Speed and 103-O) were as a rule baked for three days at 50° C., and the red-sensitive 
103-F plates were bathed in a 4 per cent ammonia solution for one minute. 


Il. IDENTIFICATION OF LINES 
Table 3 lists the lines identified with certainty in any or all of the six nebulae under 


investigation. The intensities given are relative to 78 assumed equal to 100. The meth- 
od of obtaining the intensities is described in a later paragraph. 


TABLE 2 


SPECTROGRAPH-CAMERA COMBINATIONS EMPLOYED 


Camera Focal Dispersion 

Length at Hy 
One-prism...... 3.5inches 195 A/mm 60-inch Cassegrain 
One-prism......... were 6 60 Cassegrain 
One-prism........ 10 71 60 Cassegrain 
One-prism........ Meares 18 36 60 Cassegrain 
Two-prism “VIA”: 3 200 60 Cassegrain 
Grating....... he 400 60 Cassegrain 
Two-prism “VIC”’.........| 10 66 100 Cassegrain 


All the lines identified, with the exception of those of hydrogen and helium, arise from 
forbidden transitions. The identification of the lines in Table 3 is complete with the 
single exception of \ 5670.5, which was measured on two plates of NGC 4151. For pur- 
poses of comparison, the relative intensities determined by A. B. Wyse"! for the emission 
lines in the planetary nebula NGC 7027 are given in Table 3. 

The chief features of the spectra which appear in absorption are the cyanogen band 
(effective wave length, \ 3873); H and K lines of Ca 11 (A 3934 and \ 3969); the G band, 
composed chiefly of CH (measured position, \ 4303); and the D lines of Na (X 5890 and 
5896). Since the absorption lines are usually ill-defined and sometimes absent alto- 
gether, the spectral types in Table 1 may be subject to fairly large uncertainties. 


III. RELATIVE INTENSITIES OF THE EMISSION LINES 
The conventional method of obtaining the true relative intensity of a wide emission 
line is to measure the intensity (corrected for the background) at each wave length in 
the line relative to the intensity at the corresponding wave length in a comparison star 
of known color temperature. These values are multiplied by the relative intensity of a 
black body” of temperature equal to that of the comparison star and summed over all 
wave lengths in the line. The ratio of such summations for two different emission lines, 
corrected for differential extinction, gives the ratio of relative intensities of the lines 
freed of instrumental and atmospheric absorption, i.e., the true relative intensities. 
The method described above has been used for determining true relative intensities 
for the six objects under investigation. However, the continuum of each spiral was used 


Ap. J., 95, 356, 1942. 
12 The black-body intensities may be calculated from Table VII in Jahnke and Emde, Tables of 
Functions, p. 46 (Leipzig, 1933). 
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instead of stars for comparison standards. Not only does this make it possible to use all 
the plates on which this continuum appears, but it eliminates uncertainties due to dif- 
ferential extinction, exposure time, changes during exposure of plate-sensitivity and sky 
conditions. 


TABLE 3 
INTENSITIES OF EMISSION LINES IN SIX EXTRAGALACTIC NEBULAE 


NGC 4151 | NGC 7469 | 
NGC NGC NGC | NGC | NGC 
ATOM | | | | 
1068 1275 3516 4051 | 7027* 
| Core | Wing Win | Core Wing | 
lon} 3720.77 | 80: | 140:)...... | 25:| 48: | 15: 
[Nem]......; 3968 : if 15 
He | 3970.1/ 25:7; 20: 25: 8 
Hé........| 20 10:| 25 | 20} 20/...... 35: 12 
ly... 4340.5 40 50 40t 40 35| 30] 35]...... 60:t 20 
Hem.......| 4685.8 40 25 25 | | 40 
[Arv].......| 4711.4 10 |. | 3 
H 8B. 4861.3 | 100§ 100, 100) 100 100} 100, 100 | 100 | 100 | 100 
[Om].......| 4959.5 | 400§ | 80 15) 55 | | 90 125 | 35 | 430 
(Omj.......| 5007.6 | 12008; 270) 190 | 1150]...... | 300 | 80 |1190 
[Fevu]......| 5720.9 10 | 20 | "4 
6310.2/ 85 100 30 35 15 
(O1]........| 6363.9 2 | sol...... 20 
Ha........ 6562.8» 1000§ 700 600 375 | |...... 420 
LV a] 6583.9) 50: | 200 | | 50 190 
(Sal 6734.3 140-210 40: | 180 |... 40 115 


* The intensities in the planetary nebula NGC 7027 (included for comparison) are from the paper by A. B. Wyse, Ap. J., 
95, 356, 1942. 

t The measures of \ 3869 and \ 3968 in NGC 1068 are probably affected by the strong absorptions in the underlying G-type 
spectrum. The true values are probably somewhat greater (see pp. 37, 38) 

{ The G band in the underlying continuum affects the measures of Hy in NGC 3516 and NGC 7469 (wing). The true values 
are probably somewhat larger. 

§ The effect of the narrow central minima in Ad 4861, 4959, and 5007 on the measures is negligibly small. 

|| P indicates the presence of a line which is too weak to measure. 

{ The sum of the [.V u] intensities (A 6548 and \ 6583) approximately equals that of Ha in NGC 1068. 


The method necessitates ascertaining the color temperature of the continuum of each 
nebula. For this purpose a special series of red-sensitive plates (103-F) was taken with 
the small-scale VIA spectrograph and 3-inch camera. Suitable bright comparison stars 
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of well-determined color gradient!* were taken at the beginning and end of the nebular 
exposures. The stellar exposures were lengthened to 10 or 15 minutes by interposing 
several thicknesses of copper mesh about 8 inches in front of the slit. The nebular ex- 
posures ranged from 40 to 110 minutes. A slit approximately 6 inches wide was used, and 
both nebular and stellar spectra were widened by trailing along the slit. Differential ex- 
tinction corrections were computed from Table 767 of the Smithsonian Physical Tables 
(8th ed.) by F. E. Fowle. The gradient of the nebular continuum relative to the standard 
stars was found to be closely linear from \ 4000 to \ 6700, and the black-body tempera- 
ture corresponding to this gradient was determined." 

Mean color temperatures were found in this way for the continua of NGC 1068, 3516, 
and 4151. The temperatures of NGC 1275, 4051, and 7469 were deduced by means of 
photometric comparisons with spectrograms of NGC 1068, 3516, and 4151 obtained on 
plates of the same type with similar development but taken one or more months apart. 
The adopted values of the temperatures of the continua are in Table 4. 


TABLE 4 
COLOR TEMPERATURES OF CONTINUA OF 6 EMISSION SPIRALS 
NGC T NGC T 


The investigations of Stebbins and Whitford" indicate that the continua of nebulae 
deviate from black-body radiation, particularly in the ultraviolet and infrared. I am 
very grateful to Dr. Stebbins, who, at my request, kindly measured photoelectrically the 
spectral radiation of NGC 1068 at the following mean wave lengths: AA 3510, 4220, 
4880, 5700, 7190, and 10300. He finds that the intensity distribution from 3510 to 
d 10300 is nearly identical with that of the nonemission spiral NGC 4826 (same nebu- 
lar type and spectrum as NGC 1068). From \ 4220 to A 7190 the intensity distribution 
corresponds closely to that of a black body of temperature 4800°. From \ 3510 to 
\ 4220, Stebbins’ measures indicate a temperature of 5600°. Considering the relatively 
low accuracy of my measures and the fact that the emission lines affect the photoelectric 
measures somewhat,!® the agreement with the mean value (5250°) used in this paper 
for NGC 1068 is considered satisfactory. Furthermore, even though the nebular spectral 
radiation deviates from that corresponding to the color temperatures adopted, the effect 
of such deviations or errors on the relative intensities is comparatively small. For ex- 
ample, an error of 500° in the color temperature of the continuum of NGC 1068 intro- 
duces an error at \ 6720 of 14 per cent and at A 3727 of only 9 per cent relative to the 
intensity of HB. The uncertainty is, of course, correspondingly less for lines less widely 
separated. 

The final true relative intensities are given in Table 3. The estimated uncertainty is of 
the order of 10 per cent, though it may exceed this value for NGC 7469 for which the 
measures are from a single plate. In view of the ultraviolet excess mentioned above, the 
intensities given for lines less than \ 4000 may be somewhat too small and have been 
marked with a colon. 

13 Chosen from the Greenwich list (M.N., 100, 189, 1939). The comparison stars were 13 Eri for 
NGC 1068, 21 LMi; and x Her for NGC 3516 and NGC 4151, and p; and p2 Lyr for NGC 3516. 

14 The Greenwich relative gradients must be reduced to absolute gradients by adding a correction for 
the Greenwich zero point. The value adopted is 1.10 in accordance with M.N., 100, 189, 1939. 

15 Pub. Amer. Astr., Soc., 10, 239, 1942. 


16 A rough calculation indicated that the effect of the emission lines in NGC 1068 would probably 
not be large but would tend to make the observed photoelectric color temperature too small. 
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NGC 1068.—Thirteen per cent of the total light of the nucleus in the photographic re- 
gion is in the emission lines. The intensities of the lines are similar to those in the 
planetary nebula NGC 7027 except that the [O 1], [.S m1], [N 11, and [Fe vit] lines are con- 
siderably stronger in NGC 1068. The line \ 6086 [Fe vu is unusually strong, although 
[Ca v] probably contributes to its strength. The corresponding [Ca v] line at \ 5308 is, 
however, absent or very weak. On the other hand, A 5276 of [Fe vitl is possibly weakened — 
by the superposed absorption \\ 5262-76 of Ca, Fe, and other metals. This blend would 
also probably obscure \ 5270 of [Fe 111] if it were present. Further evidence regarding 
absorptions superposed on emissions will be discussed later (p. 37). The absence of 
\ 5270 of [Fe m1] and of the lines AX 5802 and 5812 of C 1v makes it impossible to deter- 
mine with certainty whether [Fe 1] or C rv is responsible for the line \ 4658.6. 

NGC 1275.—The intensities resemble those in NGC 1068 except that the [O 1] and 
|S 11] lines are stronger and N1 and N2 [Om] considerably weaker while \ 4686 Het 
is absent altogether. 

NGC 3516.—The emission spectrum consists practically of nothing but hydrogen. 
The only other lines visible are extremely weak N1 and N2; \ 3727 [On] is absent. 
Only 5 per cent of the total photographic light of the nucleus is in the emission lines. 
The K line of (Ca 11) is visible in absorption, and the G band (A 4303) is present at the 
edge of Hy (4 4340). 

NGC 4051.—N1 and N2 are relatively weak. The plates are not sufficiently strong 
in the violet to establish the existence of \ 3727 [O m1]. 

NGC 4151.—The intensities of the cores of the lines are very similar to those in 
NGC 1068 except that the [O01] lines are considerably stronger in NGC 4151. The 
Balmer decrement of the core of the lines is the same as that of the wings. The hydro- 
gen wings contribute 8 per cent of the total photographic light of the nucleus, the rest 
of the emissions 14 per cent. 

NGC 7469.—N1 and N2 are weaker than in NGC 1068 and NGC 4151. The Balmer 
decrement seems to be normal for the cores of the lines, but that for the wings appears 
to be extremely slow. The G band and the K line of Ca are visible, despite the rela- 
tively high color temperature (8000°). 


IV. PROFILES OF THE EMISSION LINES 


All spectrograms of suitable dispersion were analyzed with the microphotometer for 
the purpose of measuring profiles of the emission lines. No attempt was made to allow 
for the instrumental profile, since only those spectrograms were used for which the 
instrumental widening was a negligible fraction of the width of the emission lines. The 
profiles of the hydrogen lines are similar when expressed in kilometers per second rather 
than angstroms. In Figure 1 are microphotometer tracings of HB, N2 and N1 showing 
the main types of emission lines observed. Figures 2-4 are the mean profiles obtained 
from the measures. Most of the lines in NGC 1068 were sufficiently symmetrical to war- 
rant combining the two halves of the line. Table 5 lists, for the lines measured in each 
object, the total widths in kilometers per second derived from the mean measured profile. 

NGC 1068.—The line widths range from 2400 km/sec for \ 4686 (He) to 3600 
km/sec for the hydrogen lines and the [.S 11] doublet \ 6717, \ 6731. The profiles of the lines 
of NGC 1068 fall into two main groups; lines with rounded maxima which are relatively 
wide in the range from 0.25 to 0.50 of the central intensity (AX 3727 [O u1],!” 4363 [O m1}, 
6086 [Fe vit],'7 6300 and 6364 [O 1], and 6717, 6731 [S m|'* and lines with sharp maxima 


17 The profile of \ 3727 is uncertain, owing to the presence of the strong absorption features at 4 3720 
and \ 3736. The profile of \ 6086 is uncertain, owing to the possibility that \ 6073 (unknown origin) and 
d 6102 [K rv] may be present. 

'S The lines \ 4069 and \ 4076 probably belong in the group, but because of blending, the profile is 
too uncertain to be included. 
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which are relatively narrow in the range from 0.25 to 0.50 of the central intensities 
(AA 3869, 3968 [Ne 111], 4686 He ur, and 4959 and 5007 [O m]). The profile of the hydro- 
gen lines seems to belong between the two groups. 

The question whether the absorption lines are also widened in NGC 1068 is not readily 
answered, because the prominent absorption lines visible (CV; K (Ca 11); Gand D) are 
either intrinsically very wide or composite. Measures from microphotometer tracings, 
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Fic. 1.—Microphotometer tracings of the emission lines \A 4860 (78), 4959 and 5007 |O 111] in the 
nebulae NGC 1068, 3516, and 4151. 


however, show that the K line in NGC 1068 is not significantly wider than the cor- 
responding line in M 32 and has about the same width as in a star of type dG5. 

NGC 1275.—Owing to the faintness of this object, only the line N1 has a well- 
determined profile. Not only is it the broadest N1 line measured, among the six nebulae 
investigated, but it shows the largest amount of asymmetry. The blue part of the line is 
60 per cent wider than the red part at J = 0.25 Imax (Fig. 4). Comparisons of the widths 
of the other emission lines from smaller-scale spectra indicate that the lines in NGC 1275 
are probably somewhat wider than the corresponding lines in NGC 1068. In particular, 
the [Ne 1] lines are unusually wide in NGC 1275. 

NGC 3516.—Only the hydrogen and N1 and N2 lines are present in emission. N1 
and N2 are weak and much narrower (1400 km/sec wide) than the corresponding lines 
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in NGC 1068 and NGC 1275. The hydrogen lines, however, are great shallow bands 
8500 km/sec wide, the widest found among the objects studied. The lines are all fairly 
symmetrical (Fig. 4). 

NGC 4051.—The N1 and N2 lines are relatively narrow (1200 km/sec), and the hydro- 
gen lines are peaked and wide (3600 km/sec). The N1, N2, and hydrogen lines appear 
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Fic. 2.—Profiles of emission lines in NGC 1068. Ordinates are intensities, abscissae distances from the 
center of each line in km/sec. The vertical lines indicate the average deviation of the individual ob- 
servations from the plotted means. For symmetrical lines the two halves of each profile have been com- 
bined. 


symmetrical within the accuracy of the measures. Small-scale spectrograms indicate that 
the remaining lines are probably similar in width to those in NGC 1068. 

NGC 4151.—The hydrogen lines are markedly different from the other lines in that 
each hydrogen line has very broad wings (7500 km/sec wide) extending outward from 
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the relatively narrow core. The discontinuity between the core and the wing is so sharp > 
that the Balmer lines appear to be composite, with sharp cores centrally superposed on 

broad emission bands. All the forbidden lines, as well as \ 4686 and the cores of the ; 
hydrogen lines, are similar in shape and relatively narrow (<1800 km/sec). ‘ 
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Fic. 3.—Profiles of emission lines in NGC 4151. Ordinates are intensities, abscissae distances from 
the center of each line in km/sec. The vertical lines indicate the average deviation of the individual ob- 
servations from the plotted means. The scale of the abscissae for the profile of hydrogen core and wing 
(lower right corner) is widely different from that of the other profiles. 


Most of the lines in NGC 4151 are asymmetrical in the sense that the blue part of the 
line is, on the average, 16 per cent wider than the red part at the point where the inten- 
sity drops to 0.25 of the central intensity. The profile of \ 6300, the only outstanding 
exception, is undoubtedly distorted by the presence of \ 6310 on the red side of \ 6300. 
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NGC 7469.—The spectrum is similar to that of NGC 4151 in that the N1 and N2 lines 
are relatively narrow (<1400 km/sec), and the hydrogen lines have cores and very wide 
wings. However, both cores and wings seem wider than those in NGC 4151 and the 
wings more intense. Absorption K (Ca 11) is extremely shallow and at least 50 A wide, 
twice as wide as in the other nebulae. The large apparent asymmetry for HB in NGC 
7469 may not be real, since the profile is determined from only one blue-sensitive plate. 
The rapidly changing background over the great width of the H@ line (8500 km/sec) 
makes the profile somewhat uncertain. 
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Fic. 4.—Profiles of emission lines in NGC 1275, 3516, 7469, and 4051. Ordinates are intensities, ab- 
scissae distances from line center in km/sec. The vertical lines indicate the average deviation of the indi- 
vidual observations from the plotted means. 


V. MINIMA IN EMISSION LINES 


An unusual feature of six of the strongest lines in NGC 1068 and one line each in 
NGC 4151 and NGC 7469 is the presence of minima in the emission lines (see Fig. 2). A 
tentative explanation of these apparent reversals is that the broad emission lines are 
superposed on a G-type spectrum which contains strong absorption lines at or near the 
position of the emission lines. Table 6 lists the emission lines in which the minima appear, 
the observed wave lengths of the minima, the suggested identifications, and the equiv- 
alent widths of the minima as measured in the nebular spectrum. Measures of the equiv- 
alent widths of the corresponding lines in the sun are included in the table for comparison. 
The measured wave lengths and equivalent widths agree with the corresponding solar 
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TABLE 5 
TOTAL WIDTHS OF EMISSION LINES (IN KM/SEC) | 
nN NGC 1068 NGC 1275 NGC 3516 NGC 4051 | NGC 4151 NGC 7469 ( 
4959, 5007. : ariel 3000 4500 1400 1200 1300 < 1400 
6300... 2700 18004 
6717\ + WOOO: 
3600 S 4500 8500 36000 7500 $500 \ 
* The true width of the line \ 6086 in NGC 1068 may be considerably smaller than that indicated, owing to the possible | 
presence of A 6073 (unknown origin) and A 6102 [K tv}. 
. dey NGC 4151, since 4 6310 influences the profile, the true width of \ 6300 is probably considerably smaller than that in- 
dicated. 
t The fact that the two lines are blended in NGC 1068 contributes only a negligible amount to the width of the line. Y 
| 
TABLE 6 
SUGGESTED IDENTIFICATIONS OF ABSORPTION MINIMA IN EMISSION LINES 
NGC Emission Observed A Suggested Identification E. W. E. W.t 
si Line of Abs. Min.* of Abs. Min. Nebula Sun 
1068. . 3869 3873 + 3872 (CN) 3 6A 4.0: A 
3968.5 (Ca 
2968 072+ ) 
1068. 3968 3972 3970 1 (He 20 10 \ 
7469. 4340.5 4300 + 4303 (CH G band) 2.0 326 
1068. 4340.5 4340 + 4340.5 (Hy) 1.0 aS 
1068. 4861.3 4862.0 4861.3 (7/73) 2.0 
| 4870.8 (Cr, Ni) 
4151. 4861 3 4873 «4871.3 (Fe) 0 6 
4872.1 (Fe) 
4957 3 (Fe) 
( 
1068... 4959 4958 1 4957 6 (Fe) 0 6 0.8 | 
5005.7 (Fe) 
5007 .3 5007 2 (Ti) 13 0.6 


5007 3 (Fe) 


* The wave length of each absorption minimum is measured on the tracings relative to the emission line 
¢t Determined from the Photometric Atlas of the Sclar Spectrum by Minnaert, Mulders, and Houtgast 
(Amsterdam, 1940). 
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values within the accuracy of the measures except for the HT (Ca 11) and He absorption 
in the \ 3968 [Ne 111] emission. The large discrepancy in this case between the observed 
and expected equivalent widths of the minimum is probably due to the great intrinsic 
width of the H/ line. The small systematic difference of approximately one angstrom be- 
tween the observed and expected wave lengths may or may not be significant. 


VI. DISCUSSION 


Ten of the twelve extragalactic nebulae known to exhibit high-excitation nuclear emis- 
sion are spirals of either early or intermediate type (i.e., Sa or Sb). The remaining two, 
NGC 1275 and NGC 3077, are objects which could be classified as either peculiar el- 
liptical or irregular nebulae. The average absolute magnitude of the twelve emission 
spirals (— 15.8) is somewhat brighter than the average nonemission spiral. All emission 
spirals have nuclei of high luminosity, which on direct photographs are scarcely dis- 
tinguishable from stars—a fact first noted by Humason in his paper on NGC 1275. 

The widths of the emission bands are probably correlated with the physical properties 
of the nucleus. The data in Table 7 suggest that the maximum width of the hydrogen 


TABLE 7* 


RELATION BETWEEN MAXIMUM WIDTHS OF BALMER LINES 
AND THE PHYSICAL PROPERTIES OF NUCLEI 


| Max. Width of Percentage of 


Abs. Mag. of | 


NGC H Lines NE ts Total Light in 
| (in Km/Sec) Nucleus 
1068... . 3600 —13.0 6 
-15.0 | 10 
7500 —14.0 | 48 
8500 —15.5 | 30: 


* The maximum widths of the hydrogen lines are from Table 5. The absolute 
magnitude and the percentage of the total light in the nucleus are derived from the 
data in Table 1. The absolute magnitude of NGC 1275 is calculated on the assump- 
tion of 05 photographic interstellar absorption, a value suggested by the low galactic 
latitude and counts of faint nebulae in the vicinity of NGC 1275. 


emission lines increases with absolute magnitude and with the ratio of light in the 
nucleus to the total light of the nebula. No definite correlation seems to exist between 
the line widths and color temperatures or spectral type of the nucleus, and the number of 
nebulae investigated is too small to indicate any certain correlations of widths of the 
emission lines with emission-line intensity ratios. 

Since some of the brightest emission knots in the periphery of other extragalactic 
nebulae have total absolute magnitudes comparable with those of the emission nuclei, the 
question arises whether such emission knots show spectra with wide lines. Spectra of 
NGC 5471 (abs. mag. —9) and NGC 604 (abs. mag. — 10), the brightest emission nebu- 
lae in the spirals M 101 and M 33, revealed no evidence of widening of the lines. At the 
writer’s request, Dr. Mayall kindly examined the available Lick Observatory spectra of 
the giant emission nebula 30 Doradus (abs. mag. — 14) in the Large Magellanic Cloud. 
He reports that the emission lines on the existing plates show no conspicuous evidences of 
broadening. Mayall also states that the very bright knot in NGC 2366 reveals no large 
line width. Evidently there is no close connection between the giant diffuse nebulae and 
the nuclei of spirals showing nebular emission lines. 
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Further discussion of the data will be presented in a contribution being prepared 
jointly by Dr. R. L. Minkowski and the present author. 


It is a pleasure to express my sincere gratitude to the National Research Council for } 
the Fellowship which made this investigation possible and to Dr. W. S. Adams, director 
of the Mount Wilson Observatory, for allowing me full use of the facilities of the observa- } 
tory. I am greatly indebted to Dr. E. P. Hubble for his continued guidance in this in- 
vestigation, and I want to thank both Dr. Hubble and Mr. M. L. Humason for generous- 
ly placing their observational material at my disposal. The helpful advice of Dr. R. L. 
Minkowski and all the other members of the Mount Wilson Observatory staff is heartily 
acknowledged. 


CARNESIE INSTITUTION OF WASHINGTON 
Mount WILSON OBSERVATORY 
October 1942 


} 

! 

} 

| 


ON THE EXISTENCE OF A THIRD COMPONENT 
IN THE SYSTEM 70 OPHIUCHI 


Dirk REUYL AND ERIK HOLMBERG 


ABSTRACT 


Photographic observations of 70 Ophiuchi at the Leander McCormick Observatory give some anom- 
alies that may be explained by the assumption of a third body. The photographic material consists 
of 33 grating plates (1932-1942) and 64 parallax plates (1914-1917, 1928-1932, 1939). The perturbations 
have a period of about 17 years and an amplitude, in X and Y, of about 07015. The mass of the invisible 
companion amounts to 0.01 solar units. 


The system 70 Ophiuchi, A and B (a = 18"0™4, 6 = +2°31’, 1900), with magnitudes 
4.3 and 6.0, is one of the most interesting and one of the most frequently observed double 
stars. The latest determination of the orbit by K. A. Strand? resulted in a period of 
87.85 years and a semi-major axis of 47558. Certain apparently periodic deviations in 
the older visual observations have at various times given rise to the suspicion that the 
system contains a third, invisible body, the period of the perturbation ranging from 6 to 
36 years. However, these observations have a comparatively low accuracy, and it is 
probable that the observed deviations from Keplerian motion must be attributed to the 
effect of systematic errors in the visual observations rather than to the influence of a 
third component. In the investigation by Strand, photographic observational material 
of high accuracy is used for the first time. According to the author, these observations 
do not support the theory of a third body. However, the photographic observations 
discussed by him extend over comparatively short intervals of time, 1915-1922 and 
1931-1935. A discussion of all the material up to and including the 1941 McCormick 
observations was presented at the sixty-eighth meeting of the American Astronomical 
Society at the Yale Observatory in June, 1942. In the present paper a new discussion 
will be made, based on Strand’s material together with the photographic observations 
made at the McCormick Observatory during the period from 1914 to 1942. The former 
consists of 29 grating plates taken at the Potsdam and Johannesburg observatories, the 
latter comprises a total of 1150 measured exposures on 97 grating and parallax plates. 

On account of their relatively low accuracy, the older visual observations will not be 
taken into account. 

The McCormick observations of 70 Ophiuchi are listed in Table 1. Most of the mate- 
rial consists of objective grating plates taken by one of us during the period 1932-1942. 
A minor part is represented by 64 parallax plates secured in the course of the systematic 
parallax work at the McCormick Observatory. For the grating plates three different 
kinds of gratings have been used, giving reductions of 2™9, 2™5, and 1™8, respectively. In 
the last case the secondary image of the bright component is practically identical in 
magnitude with the primary image of the faint component. Thus, possible guiding errors 
are effectively neutralized. The grating plates have been taken only under the best pos- 
sible conditions of seeing, with the result that the quality of nearly all the images is good. 
The exposures, ranging in number from 15 to 60 per plate, are placed in a regular pattern, 
close to and symmetrically located with respect to the center of the plate. Only the best 
images were selected for measurement. The number of measured images, which averages 
about 30 per plate, is given in the second column of the table. The measurements have 
been made by the present writers, each of them accounting for about half the number of 
plates. In order to check possible systematic differences, we have measured a number of 
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Epoch 


1917.47... 
1929.54... 
1931.70. 
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TABLE 1 


42°78 
+3.436 
45.454 


XQ 


oN 


6 066 
+6 .067 


6.070 
6 056 
6.073 
6.073 
+6 .069 


ICCORMICK MATERIAL 


+0006 
006 
007 
007 


.006 
003 


.010 
004 


005 
005 
005 
003 


007 
004 
012 
O15 
003 


005 
010 
006 

OOF 


006 
O11 
007 
O04 


005 
005 
O04 


005 


004 
006 
005 
O06 
003 


003 
005 
005 
0.002 


| —3"469 +.0°005 
| 3.677.005 
—3.592 .005 
005 
3.367 006 
3.447 004 
3.431 .006 
3.429 004 
— 3.424 002 
3.301 .004 
3.318 009 
3.292 O04 
—3.302 003 
3.222 
3.205 005 
3.207 005 
3.218 004 
2008 
3.422 004 
3.124 004 
3.130 013 
3.135 005 
3.146 009 
— 3.129 O02 
3.057 005 
3 060 008 
3.042 005 
2 963 007 
2.969 O11 
2.941 007 
2.951 006 
| —2 956 
2.842 004 
2.824 005 

| 

| | 

| —~2816 | 
2 650 .004 

| 2.636. | 

2 635 O06 | 

| 2.630 .006 | 

| —2 639 003 | 

| | 

2 504 007 | 
2.491 005 
2 508 004 

—2.512 0.004 | 


Remarks 


11 paral. pl. 
9 paral. pl. 
13 paral. pl. 


13 paral. pl. 


No grat. 
Grat. 2.9 
No grat. 
Grat. 2.9 


Grat. 2.5 
Grat 2.5 
Grat. 2.5 
Grat. 2.5 
Grat. 2.5 
Grat. 2.5 
Grat. 2.5 
Grat. 

Grat. 

Grat. 


Grat. 
Grat. 


Grat. 1.8 
Grat. 1.8 
Grat. 1.8 


Grat. 1 
Grat. 1 
Grat. 1 
Grat. 1 


Grat. 1 
Grat. 1 


18 paral. pl. 


Grat. 1 
Grat. 1 
Grat. 1.8 
Grat. 1 


Grat. 1 
Grat. 1.8 
Grat. 1 
Grat. 1.8 


— 


) 
| No. X p—a (2000.0) 
18 
25 | 
| | 
30 / 
33 | 
25 
33 
+ 
32 ) 
43 
22 
| 
42 
12 
1936.59... 
15 6 024 
See 40 6 038 
+6 040 
eer, 22 6.053 8 
19 6 067 8 
33 6 055 8 
19039200 | 27 6 080 | 
aural 27 6 082 
+6 081 
| 
1939.76... | 36 +6077 
6.054 
1942.47.........| 53 
40 
| 
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plates in common. However, no differences of a systematic nature could be found. The 
plates have been measured in the usual way on the two Gaertner measuring machines of 
the McCormick Observatory. Each plate has been oriented by means of a trail of a 
near-by bright star and measured in direct and reversed positions. 

In Table 2 the deviations of the observed values from those corresponding to the 
latest orbit, as derived by Strand, are listed. The positions corresponding to the orbit 
have been derived graphically, which gives an accuracy suflicient for the present purpose. 
In order to have a complete picture, the seven normal places corresponding to Strand’s 
observations have been included in the table. A superficial inspection of the table shows. 
that the deviations are of a systematic nature. This is especially the case in declination 


TABLE 2 
DEVIATIONS OF OBSERVED VALUES FROM STRAND’S ORBIT 


Ricut ASCENSION DECLINATION 
Epocu Wr. MATERIAL 
Obs. Dev. Obs. Dev. 
1914/80) 1 +2"778 —("029 —3"469 | +0028 | McCormick paral. pl. 
15.25. 3 2.915 | — .005 3.521 | + 011 | Strand 
16.80..... 3 3.298 | + .005 3.637 — .004 Strand 
Vc. eo 1 3.436 | — .010 3.677 — .010 | McCormick paral. pl. 
1925).3 4 3.837 + .002 3.744 | — .008 Strand 
I 2 4.438 | — .O11 3.802 | — .017 Strand 
29.54... 1 5.454 | — .007 3.592. | + .012 | McCormick paral. pl. 
31.27 7 5.037 | 3.497 | + .011 | Strand 
31.70... 1 5.674 | — .002 3.482 | — .001 | McCormick paral. pl. 
a 6 5.734 | — .007 3.428 | + .005 | Strand 
32.47. ye 2 5.711 | — .031 3.424 | + .009 Reuyl 
2 5.906 + .019 3.302 | — .016 | Reuyl 
S556: «.; 5 5.960 + .010 3.202 — .003 | Strand 
5.939 — .O13 | — .016 | Reuyl 
36.59. 3 5.990 — .009 3.129 | — .018 | Reuyl 
37.58 3 6 040 + .003 3.052 | — .032 | Reuyl 
38.56 3 6 061 — .006 2.956 | — .028 | Reuyl 
39.60 3 6 081 — 007 2.833 — 010 | Reuyl 
39.76 1 6.077 — O15 2.816 — .010 | McCormick paral. pl. 
41.59. 3 6.067 — .036 2.639 — .027 | Reuyl 
1942.59 3 +6 .069 —0.025 —0.013 Reuyl 


where the last ten deviations all have a negative sign, their average value amounting to 
about —0"017. For the subsequent least-squares solutions the weights given in column 
2 of the table have been adopted. In order to take into account the different gratings 
used and, furthermore, to avoid an excessively large range in weights, a simple system of 
weighting, equal to w/P.E., has been adopted. In the case of plates taken without grat- 
ing the quantity w is taken as 07005, whereas for grating plates w is put equal to 07007, 
0” 008, and 0010, corresponding to gratings of 2™9, 2™5, and 1™8 reduction, respectively. 
In this way the normal places obtained from parallax plates are given a weight of 1, 
whereas most of the normal places corresponding to the McCormick grating plates re- 
ceive a weight of 3. 

Since the present material is still not large enough to justify a new determination of 
the orbit of components .1 and B, the deviations listed in Table 2 will be used as a start- 
ing-point for the study of possible perturbations. For a least-squares solution, the follow- 
ing relation will be used: 


A=a sin 
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The expression a@ sin 6 (¢ — fo) corresponds to a motion in a circular orbit by the third 
body. The present material obviously is not large enough to reveal a possible eccen- 
tricity. The amplitude of the perturbation equals a, whereas 0 is inversely proportional 
to the period (6 = 27/P). The term c + (¢ — ¢éo) » gives a correction to Strand’s orbit. 
To make possible a least-squares solution the above expression is differentiated: 


dA = 


where 
A=sin b(t-—t), B=al(t—t) cosb(t—t), 


Preliminary values of a, 5, and ¢9 are derived graphically in order to determine the fac- 
tors A, B, and C. A least-squares solution then gives the corrections da, db, dto, dc, and 


du. 


25 35 40 
| 
ee 
e e 


+'03 4 


1915 30 40 


Fic. 1.—Deviations of observed values from Strand’s orbit 


By using the above procedure we obtain the following results: 


RIGHT ASCENSION DECLINATION 
a= 0014 + 07003 a= 0"015 + 07002 
P= 1673 0x9 P= 1779 
to = 1930.8 to = 1925.6 
C= —O0*011 0”002 C= 0”003 
0”0000 070002 —0"0005 070002 


Thus, a perturbation is obtained with a period of about 17 years and an amplitude, in V 
and Y, of about 07015. The residuals obtained in the least-squares solutions correspond 
to probable errors of 070075 in right ascension and 070051 in declination, the average 
weight being 2.86. It may be remarked that the average observed probable errors of the 
positions of Table 2 amount to 070031 and 070028, respectively. 

Figure 1 gives an illustration of the anomalies found above. By reference to Table 2 
it will be seen that the normal points of greatest weight are represented in Figure 1 by 


70 OPHIUCHI 45 


the filled circles of largest diameter. The perturbation seems to be rather well established 
in declination, whereas in right ascension additional observations are desirable. In the lat- 
ter case the form of the curve depends to a large extent on the observations made during 
1941 and 1942. Additional observations during the next few years will decide whether 
or not the right-ascension curve has to be changed. 

Adopting a parallax of 0720, the above amplitude of 07015 corresponds to 0.075 
astronomical units. The number of parallax plates available at the present time is in- 
sufficient to determine whether the third body belongs to component A or to component 
B. In the first case, if the mass of component A is assumed to be equal to 1.1 solar units,? 
a period of 17 years would correspond to a total orbital radius of 6.8 astronomical units. 
Thus, the resulting mass of the third component would be 0.012 solar units. In the second 
case, if the mass of component B is taken as 0.7 solar units, the total orbital radius would 
be 5.9 astronomical units, resulting in a value of 0.008 solar units for the mass of the 
third body. 

As to the possibility of the blending of a fairly bright third component with A or B, 
this has been dismissed in view of the value of the orbital radius of 174 or 172. It may be 
remarked that a spectroscopic study of a perturbation of this kind would be extremely 
difficult, since the orbital velocity of the star (component A or B) would amount to only 
about 0.15 km/sec. The orbital radius of 6 or 7 astronomical units is rather large, as com- 
pared with the semi-major axis of the orbit of components A and B of about 23 astronom- 
ical units. It is to be expected that the orbit of the third body will not be very stable. 
The period, as well as the amplitude, of the observed perturbation may change consider- 
ably in the course of time. When more extensive photographic observational data are 
collected in the future, the system 70 Ophiuchi may furnish an illustration of a very 
interesting three-body problem. 

LEANDER McCorMICK OBSERVATORY 


UNIVERSITY OF VIRGINIA 
October 1942 


2 Holmberg, Determination of the Mass Ratio of 70 Ophiuchi (in press). 


3 The following photographic observations of high quality, listed in Pub. A.A.S., 10, 70, 1940, have 
been inadvertently omitted from the discussion and should be entered in Table 2: 


Ricut ASCENSION | DECLINATION 
Epocu | MATERIAL 
Obs. Dev. | Obs. | Dev. 
1937.62.........| +67026 —0'012 —3*020 —0'004 Lick 
| +6.084 — .601 —2.851 — .017 Sproul 
| +6.078 —0.023 —2.732 —0.010 Sproul 
| 


The average deviation of these values from the curves is 0"000 in right ascension and +0"012 in 
declination. Their inclusion will result in a small change in the curve in declination but will leave the 
curve in right ascension unchanged. It is noted that there will be no effect upon the discussion of the 
masses in the system. 
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A SPECTROGRAPHIC STUDY OF METEORITES 
W. W. A. JoHNSON AND DANIEL P. NORMAN 


ABSTRACT 

Eighteen meteorites were investigated spectrographically for the presence of sixty-nine elements. 
Forty elements were found, twelve of which had not previously been reported by King. Thespectrograms 
were taken by a jumping-plate technique that permitted estimates to be made of the relative volatility 
of the elements found, all of which, with the exception of aluminum and titanium, vaporized in the order 
expected from considerations of their boiling-points and abundances. No obvious grouping of these 
meteorites could be made on the basis of their atomic constituents. The two iron meteorites investi- 
gated contained less titanium than the stones, in agreement with previous chemical and spectroscopic 
data. 


Dr. A. S. King’s! interesting spectroscopic examination of thirteen meteorites sug- 
gested a further spectrographic study stressing the differences rather than the similarities 
in the minor atomic constituents of meteorites. Mr. H. H. Nininger, of the American 
Meteorite Laboratory, very kindly furnished us with samples of eighteen meteorites, 
which he described as follows (they have been arbitrarily numbered by us): 

. Holbrook, Arizona; fell July 19, 1912 (described by Foote and Merrill, Smithsonian Misc. 
Coll., 69, No. 9, Pub. 2149) 

2. Miami, Texas; date of fall unknown 

3. Washougal, Washington; July 2, 1939 (sawings produced by sawing very friable material) 

4. Whittman, Nebraska; date of fall unknown 

5. Johnstown, Colorado (crystalline inclusion separated from matrix) 

6 

7 

8 


— 


. Springwater, Saskatchewan; date of fall unknown (pallasite) 

. Loomis, Nebraska; date of fall unknown 

. Texline, Texas; date of fall unknown 
9. Arcadia, Nebraska; date of fall unknown 
10. Smith Center, Kansas; date of fall unknown 
11. Oxidized crust of Monahans, Texas, iron meteorite 
12. Richardton, North Dakota; fell June 30, 1918 (Jour. Geol., 27, 431) 
13. Sandia Mountains, New Mexico; date of fall unknown 
14. Pasamonte, New Mexico; fell March 24, 1933 
15. Leeds, Quebec (nodule) 
16. Johnstown, Colorado; fell July 6, 1924 (described by Hovey and Merrill) 
17. Sofer, Oklahoma; sawings from an iron meteorite 
18. Graphite nodule from Xiquipilco, Mexico, iron meteorite 

(Meteorites Nos. 12, 14, and 18 appear to be common to King’s list and ours) 


The spectrographic technique used has been developed by the authors on the basis of 
qualitative analyses of a wide variety of ferrous samples of known composition. A pair 
of special spectrographic electrodes of quarter-inch diameter (Dow Chemical Co., Mid- 
land, Michigan) were introduced in the arc stand. In the lower (cathode) electrode a 
crater 5 mm wide and 6 mm deep was drilled by a special fixture. The upper electrode 
was pointed in a pencil sharpener reserved exclusively for this purpose. The arc was 
operated at 15 amperes (d.c.) for 1 minute to volatilize any traces of metal that might 
have been introduced in preparing the electrodes; then a 40-second spectrogram was 
taken to identify any impurities intrinsically present in the electrodes (boron was usually 
the only impurity; occasionally there were traces of magnesium and silicon), and the elec- 
trodes were allowed to cool in the arc stand. The sample (weighing approximately 100 
mg.) was introduced in the crater of the lower electrode, the spectrograph shutter opened, 
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the arc started, and the spectrograph plateholder racked down every 40 seconds until 
the sample was burned to completion. (This procedure is essentially the jumping-plate 
system used in eclipse spectrography.) The exposure time used was the optimum ex- 
posure dictated for our instruments by the twin opposing factors of maximum exposure 
for faint lines and minimum exposure for the background to permit maximum visibility 
for faint lines. 

The direct-current arc was operated at 15 amperes (230-volt input), while the inter- 
electrode gap was held at 3 mm as the carbons burned away. The gap length and am- 
perage were maintained rigidly. The voltage across the gap fluctuated roughly from 40 
volts at the start of a run to 50 volts when the sample had been completely consumed. 
Even the most refractory samples were completely volatilized at the end of 10 minutes’ 
arcing. We cannot stress too strongly the importance of thus burning samples to comple- 
tion, since small traces of refractory elements (particularly the noble metals) and of ele- 
ments which form high-boiling carbides in the reducing atmosphere present in the arc, 
notably zirconium,? will be markedly concentrated in the last few minutes of exposure 
time. 

An enlarged image of the arc was formed on the slit of the spectrograph with a quartz 
lens and was so adjusted that the image of the cathode fell just off the edge of the Hart- 
mann slide delimiting the slit length, thus making use of the enhanced sensitivity of 
many elements in the cathode layer.’ Two grating spectrographs were used in making 
the spectrograms. The wave-length range from \ 2175 to \ 3525 was photographed on 
Eastman 33 plates in the first order of a 3-meter grating spectrograph having a disper- 
sion of 5.6 A/mm. The wave-length region from \ 3200 to \ 6800 was photographed on 
Agfa Superpan Press film in the first order of a 1-meter grating spectrograph having a dis- 
persion of 15A/mm (the second-order ultraviolet was absorbed by a Wratten 2A filter). 

The following sixty-nine elements (which are usually detected in qualitative spectro- 
chemical analyses) were sought: Ag, Al, As, Au, B, Ba, Be, Bi, Ca, Cb, Cd, Ce, Co, Cr, 
Cs, Cu, Dy, Er, Eu, Fe, Ga, Gd, Ge, Hf, Hg, Ho, In, Ir, K, La, Li, Lu, Mg, Mn, Mo, Na, 
Nd, Ni, Os, P, Pb, Pd, Pr, Pt, Ra, Rb, Re, Rh, Ru, Sb, Sc, Si, Sm, Sn, Sr, Ta, Tb, Te, 
Th, Ti, Tl, Tm, U, V, W, Yb, Yt, Zn, Zr. The elements found are listed in Table 1. The 
first column of the table contains the elements in the alphabetical order of their symbols. 
The succeeding eighteen columns list the approximate abundance of the elements found, 
on the following logarithmic scale: 


1 = 100-1 per cent 5 = Less than 0.01 per cent 

2= 10-0.1 per cent ? = Abundance estimate uncertain 

3=  1-0.01 per cent m = Partially masked by lines of other constituents. 
4=  0.1-0.001 per cent 


The last two columns list the elements found by King,‘ in irons and in stones, respective- 
ly. His numerical intensity scale has, however, been inverted to agree with our magni- 
tude scale. 

The limit of identification® for the spectrographic analysis of the elements varies mark- 
edly with the element concerned and with the matrix in which it is present.® Since iron is 
a major constituent of all meteorites, stones as well as irons (the fact that carbon is the 
major constituent of graphitic meteorites is counterbalanced by the fact that all the 
samples are burned in graphite electrodes), the matrix may be considered, to a first 


2D. Richardson, Proc. Fifth Summer Conference on Spectroscopy, p. 68, New York: John Wiley & 
Sons, 1938. 


3L. W. Strock, Spectrum Analysis with the Carbon Cathode Layer, London: Adam Hilger, 1936. 
‘Op. cit., p. 510. 

5F. Feigl, Mikrochem., 1, 4, 1923. 

® Pierce, Torres, and Marshall, Ind. Eng. Chem. Anal. Ed., 12, 44, 1940. 
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TABLE 1 
ABUNDANCES OF ELEMENTS IN METEORITES 
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approximation, as reasonably similar in all the meteorites and as comparable to terrestrial 
ferrous samples. However, even in the case of samples having identical major constit- 
uents, minor constituents, such as sodium, can be detected in concentrations of 1x10~° 
per cent, whereas an element such as tungsten might not be detected in concentrations 
of 1X10~? per cent. 

Our estimates of the abundances of the various elements found have been corrected 
for these differences in sensitivity on the basis of our experience with samples of known 
composition. However, such a correction cannot be applied to the twenty-nine elements 
(out of the sixty-nine sought) which were not detected. These elements can be divided 
roughly into two groups: elements whose limit of detection is of the order of 10~* per 
cent or greater and elements whose limit of detection is 10~* per cent or less. These two 
groups are listed respectively in the first and second columns of Table 2. The third col- 
umn lists the elements not found whose ultimate lines may have been masked by sub- 
ordinate lines of other elements. 

TABLE 2 


SENSITIVITY OF ELEMENTS NOT FOUND 
IN METEORITES 


Less than Masked 
107" Per 10° * Per Cent 
Bi Dy Nd Tb B 
Cb Er Os Th Ga 
Ge Eu Pr Tl Gd 
Ra Hf Re U Rb 
Rh Ho Sm W 
Vd La Ta Yt 
Lu 


The following twelve elements not reported by King‘ were present in the meteorites 
we analyzed: As, Be, Cd, Hg, In, Ir, Pd, Pt, Te, Tm, Zr. We note that they are all 
present in amounts of four or less. King found three elements in the meteorites he studied 
that we did not find: Rb, Ga, Ge. The ultimate lines of rubidium (A 7800 and \ 7947), 
used by King, were not on our plates, and the lines we used (A 4202 and 4215) are 
masked when iron is a major constituent. The ultimate lines of gallium were masked by 
subordinate lines of other elements. The ultimate lines of germanium were favorably 
placed in the wave-length range covered by our plates and would not have been masked 
by lines of other elements. Repeated search failed to disclose any trace of germanium, 
and we must conclude, therefore, that it simply was not present in the samples we 
analyzed. 

The various elements observed should appear and disappear, in successive spectra 
taken by the jumping-plate technique, in the order of their boiling-points and abun- 
dances. Measurements show that the elements found fall into four groups: (1) elements 
which appear only in the first two spectra: Ag, Ba, Ca, Cd, Cs, Cu, Hg, In, K, Li, Na, Pb, 
Sn, Sr, Zn; (2) elements persisting about halfway down the plates: Be, Mg, Mn, Mo, P, 
Si; (3) elements appearing in the last few spectra only: As, Au, Ir, Pd, Pt, Ru, Sb, Te, 
Tm, Zr; and (4) elements persisting throughout all spectra: Al, Ce, Co, Cr, Fe, Ni, Ti, V. 
We note that aluminum and titanium in terrestrial samples usually occur in group 2; and 
an explanation of their anomalous behavior would probably have to be sought in the 
nature of the negative radicals attached to them, since the latter are known to affect 
their volatility.7 Elemental phosphorous boils at about 270° C.; but anions containing 


7 Brode and Hodge, J.0.S.A., 31, 58, 1941. 
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phosphorous dissociate at very much higher temperatures, even in a reducing atmos- 
phere, and hence the appearance of phosphorous in group 2 is not unexpected. In gen- 
eral the time sequence of the elements is in good agreement with expectations based on 
their volatilities and abundances. 

An attempt was made to group these meteorites on a rational system based on the 
observed abundances of various elements. We note that iron was always high and could 
not be used for classification, in agreement with King’s observation that iron was high in 
stones, irons, and graphite nodules. 

Titanium, however, appeared to offer a satisfactory peg on which to hang a classifica- 
tion. Accordingly, the meteorites were divided into two groups depending on (1) 
whether titanium was abundant or of intermediate abundance or (2) whether titanium 
was poor or lacking. An attempt to divide (1) further on the basis of aluminum content, 
yielded two subdivisions: (A) high in A/ and 77: Nos. 3, 4, 5, 9, 14; (B) low AI, inter- 
mediate 7i: Nos. 1, 2, 7, 10, 12, 16. A different regrouping of these meteorites rich in 
titanium was attempted on the basis of the presence or absence of certain volatile ele- 
ments (As, Li, Sr, Zn), which yielded (A’) volatiles present: Nos. 4, 5, 10, 12, 14; and 
(B’) volatiles absent: Nos. 1, 2, 3, 7,9, 16. We do not feel that great significance should 
be attached to these correlations, and no further correlation is apparent in the data for 
this group. 

The meteorites low or lacking in titanium (2) could in turn be subdivided into two 
groups depending on their silicon content: (C) high silican: Nos. 6, 8, 15, and (D) low 
silicon: Nos. 11, 13, 17, 18. We note further that in these groups high silicon goes hand 
in hand with high magnesium, while low silicon is accompanied by intermediate mag- 
nesium content. There is a virtual absence of aluminum and sodium in group D, while 
phosphorous is higher in group D than in group C. 

Considering the diversity of origin and of physical structure of the samples we have 
studied, it is not surprising that they fail to fall into simplified groupings. 

An important question remains unanswered at present: What are the variations in the 
distribution of the minor elements throughout a single meteorite, and is there any evi- 
dence of a radial inhomogeneity in the distribution of the more volatile elements as we 
go from the center of a meteorite outward? We know that the inclusions in meteorites 
differ markedly from the matrix (cf. Nos. 5 and 16), and we know that, mineralogically 
speaking, meteorites are not uniform (pallasite, e.g., consists of a magnesium-iron silicate 
in a network of metal*); but we know of no studies of the variation of representative 
samples taken along radial borings of a single meteorite. We hope, in the future, to be 
able to carry out such studies. 
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THE CONTINUOUS SPECTRUM OF MODEL 
STELLAR ATMOSPHERES 


Racpu E. WILLIAMSON 


ABSTRACT 


On the assumption that atomic hydrogen and the negative hydrogen ion are the only sources of con- 
tinuous opacity in stellar atmospheres, the properties of the continuous spectra of twelve model atmos- 
pheres characterized by various values of the effective temperature 7, and the average electron pressure 
pe, are investigated. An improved value of the atomic continuous absorption coefficient for H~ is used in 
the computations. Approximate values for the Rosseland mean absorption coefficients « are derived for a 
wide range of the parameters p, and 7, and for the twelve model atmospheres (with 10 < p, < 104; 5600° 
< T, < 10,080°); the monochromatic mass-absorption coefficients x, are calculated for a number of se- 
lected frequencies covering a wave-length range of 2500—10,000 A. The ratio x/xy is used to determine the 
logarithm of the ratio of the monochromatic flux /, through each model atmosphere to the black-body 
flux © corresponding to 7, for the atmosphere. Color temperatures for regions 3300-3650 A, 3850- 
4550 A, and 5550-7150 A, as well as the magnitudes of the Balmer and the Paschen discontinuities in the 
continuous spectrum, are computed for each of the model atmospheres. Comparison of the calculated val- 
ues of the Balmer discontinuity with observations of this quantity indicate a value of p, = 10° for stars of 
spectral types A and a somewhat smaller value for stars of type F. Plots of the reciprocal color tempera- 
ture @. = 5040/7, in three spectral regions, as a function of @, = 5040/7, for observations on stars, for 
the models here considered, and for those of Wildt, who used older values of the absorption coefficient of 
H~ (both the latter for p, = 10%), show that the agreement of theory and observation has been improved 
by the present calculations. In general, the agreement with the observations is less satisfactory in the 
longer wave lengths. It is shown that the present models predict some increase of T. with frequency for the 
visual and near infrared, while the observations of Williams and of Hall show a slight decrease of T, with 
frequency. The magnitude of the Paschen discontinuity for Vega as observed by Hall and Williams is 
satisfactorily predicted by an appropriate model (0, = 0.5, p. = 10%). 


I. INTRODUCTION 


An important problem of theoretical astrophysics is the prediction from purely 
physical principles of the radiant flux F, emerging from the surface of a star. While an 
ideal treatment of the problem would involve an accurate computation of F, (including 
the rapid variations in it which will constitute the absorption lines) for a number of model 
atmospheres defined by certain values of the parameters 7., g, and the chemical abun- 
dances, it appears that in practice we can formally distinguish between a problem of the 
continuous spectrum of a star and a problem of the absorption lines. 

On this understanding we suppose that we can consider the continuous spectrum of a 
star as arising from sources of opacity which change only slowly with frequency and can 
treat absorption lines as superposed on this background in such a way that the original 
large-scale distribution of energy is unaffected by them. 

The calculation of the energy distribution in the continuous spectrum of a model 
stellar atmosphere of specified constitution was first undertaken by McCrea.' He investi- 
gated atmospheres composed of a mixture in different proportions of hydrogen and a 
hydrogen-like element with an ionization potential of 5 volts and compared the calcu- 
lated continuous spectra with the observations of Yii. Since then similar calculations 
have been made, assuming various chemical compositions of the model atmospheres, by 
Unséld? and by Pannekoek.* However, Wildt’s discovery‘ of the importance of the nega- 
tive hydrogen ion as a source of opacity in stellar atmospheres has necessitated a revision 


1 M.N., 91, 836,1931. 
2 Physik der Sternatmos pharen, Berlin: Julius Springer, 1938. 
% Pub. Astronomical Institute of U. of Amsterdam, No. 4. *ApJ., 89, 295, 1939. 


51 


1OS- 
en- 

on 

ld 
uld 
in 

Ca- | 
(1) 
im 
nt, 
er- | 
in | 
le- | 
nd 
Id 
or 
vO 
id 
a 
1e 
ve 
1e 
4S 
y 
e 
e | 
e | 

|_| 


52 RALPH E. WILLIAMSON 


of these calculations. For the continuous spectrum this has been done by Wildt,® while 
Strémgren® has incorporated it into his treatment of the line spectrum of the sun. 
Wildt’s work was based on values of the continuous absorption coefficient for H~ as com- 
puted by Massey and Bates,’ who used an approximate wave-function for the ion derived 
by Bethe® and by Hylleraas.* It has been recently found by the writer!® that a more 
accurate wave-function for H~ yields an absorption coefficient substantially different 
from that obtained by Massey and Bates. It is therefore of some interest to see in what 
way the change in the absorption coefficient will affect the calculations of Wildt on the 
continuous spectra of model stellar atmospheres. This is the principal object of this 
paper. More particularly, we shall investigate the emergent flux F, of the radiation in 
the continuous spectra of model stellar atmospheres, whose continuous opacity is sup- 
plied entirely by atomic hydrogen and negative hydrogen ions."' In the calculations the 
customary assumptions of local thermodynamic equilibrium and of the constancy of the 
ratio x,/x throughout the atmosphere have been made. Because of the latter assumption 
it is convenient to use as the independent parameters the effective temperature 7", and 
the electron pressure p.'? (instead of g). 


II. THE ROSSELAND MEAN ABSORPTION COEFFICIENT 


It is well known from the classical investigations on radiative equilibrium of stellar 
atmospheres by Milne and others'* that the emergent flux F, is given by 


Ky (+-) 
(1) 


C2 — 1 x 


where x, and x denote the monochromatic and the Rosseland mean mass-absorption 
coefficients, respectively. Further, in equation (1) 7 denotes the optical depth derived 
from x and g(r) is a slowly varying function of + and approximately equal to 3. Finally, 


= (2) 
1 


Burkhardt" has evaluated equation (1) numerically and gives a table of the quantity 
logio F,/F‘® as a function of the two parameters hv/x7,. and «/x,, where 


(3) 


5 ApJ., 90, 611, 1939; 93, 47, 1941. 
5 Festschrift fiir Elis Strémgren, p. 218, Copenhagen, 1940. 
7Ap.J., 91, 202, 1940. ° Zs. f. Phys., 60, 624, 1930. 
§ Zs. f. Phys., 57, 815, 1929. 10 Ap J., 96, 438, 1942. 
' For a discussion of the validity of this assumption cf. Wildt, Ap./., 90, 611, 1939. 


'2 Thus the presence of free electrons from other sources than hydrogen is assumed, but no calculations 
based on abundances and ionizations of the ‘“‘metals’’ are considered. Such calculations more properly 
belong in the theory of the formation of absorption lines (see Strémgren). 


13 Cf. Handbuch der Astrophysik, 3, Part III, 130, 147-55, Berlin: Julius Springer, 1930. 
14 Zs. f. Ap., 13, 36, 1936. 
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It is therefore necessary to have a knowledge of both x, and x for each pair of values of 
T, and p. for which we desire the emergent flux /,. Since the run of F, with frequency is 
insensitive to small errors in « and since a direct evaluation of k would be quite tedious, 
the following approximate method may provide sufficient accuracy for our present pur- 

oses. 

. Let the Rosseland mean of the mass-absorption coefficient of atomic hydrogen alone 
be x(H) and that due to the negative hydrogen ion alone, x(H~). For most values of 7, 
and p, one of the two mean absorption coefficients is much larger than the other. We 
can then take account of the presence of the less important source of opacity as a small 
correction to the more predominant one. Such a table of corrections has been pro- 
vided _by Strémgren" in his work on the solar atmosphere. He tabulates the quantity 
logio k/K(H) as a function of logio K(H- )/x(H) and of 7. It is further seen that the de- 
pendence on 7 of this correction term is practically negligible. The problem of obtaining 
the Rosseland mean x is now greatly simplified, involving only the calculation of x(H~) 
and applying the appropriate correction term to the already well-known quantity x(H).!® 

The contribution to the monochromatic absorption coefficient of H~ is made by 
bound-free and free-free transitions of the optical electron. The latter is comparatively 
unimportant except for wave lengths greater than 12,000 A. 

Let a,(H—) denote the absorption coefficient per negative hydrogen ion due to the 
bound-free transitions; further, let a{/)(H~) denote the absorption coefficient per neutral 
hydrogen atom and per unit electron pressure due to the free-free transitions. Then 


x» (H-) = at — | 


my 


Nu- = Nu 


a.) (H-) 


N, denotes the number of X-atoms (or X-ions) per unit volume, and my is the mass 
in grams of the hydrogen atom. Now 


Nube _ 4 (2am) (kT) = , 


(5) 
3/2 5/2 


where hyx is the first ionization energy of the X-atom or X-ion. Since Ny-< < Nx for 
all cases of interest in stellar atmospheres, we may write equation (4) as 


The effect of stimulated emissions (as negative absorptions) is readily taken into ac- 
count. Denoting the “effective” absorption coefficient by x,, we have 


= Ky 


(7) 


15 Op. cit., Table 4. 
'6 See, e.g., Wildt, Ap.J., 90, 611, Table 1, 1939. 


| 
| 
4 
y 


54 RALPH E. WILLIAMSON 


Now 
a(H-) (8) 
Mr bot pe 
where 


In future by “absorption coefficient” we shall always mean the “effective”? absorption 
coefficient including the stimulated emissions as in equation (7); we shall accordingly 
suppress the primes. Since the integrand in equation (9) is independent of p,, it is clear 
that one set of integrations for different values of 7 is all that is needed to obtain 
x(H~-). Table 1 gives d(H~) as computed by graphical integration, using the revised 


TABLE 1 
ROSSELAND MEAN OF H~ ABSORPTION COEFFICIENT/GRAM NEUTRAL H/DYNE /, 


Wheeler and Wheeler and 
P | eele 
i & resent Wildt 7 Present | Wildt 
14.5 12.4 {| 3.5 
21.8 19.4 23.1 23.9 


values of Williamson!’ for a,(H~) and the results of Wheeler and Wildt for a‘! (/7~).!8 
The last column of Table 1 gives for comparison the results of a similar computation by 
Wheeler and Wildt, who used the earlier values of the bound-free absorption coefficient 
according to Massey and Bates. It is seen that no large changes in d(//~) have been in- 
troduced by the new calculations. By interpolating among the values given there and 
substituting into equation (8), we can readily obtain k(H~). In Table 2 the common loga- 
rithm of this quantity has been entered in the top line of each square; further in this 


table the second line contains logig x(/7) from Wildt’s paper,'* while the third gives logy, k, 
derived after applying Strémgren’s corrections. 


Ill. THE CONTINUOUS SPECTRUM OF THE MODEL ATMOSPHERES 


As we have already stated in the introductory section, to determine the emergent flux 
F, for a given model stellar atmosphere, we need the value of the absorption coefficient at 
the particular frequency v in which we are interested, as well as the Rosseland mean 
absorption coefficient x. We have already provided the necessary table of « in the pre- 


17 Loc. cit. 
18 281, 1942. 
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ceding section. It now remains to evaluate x,, which is the sum of the absorption coef- 
ficient x,(H~) due to H~ and the absorption coefficient x,(H) due to hydrogen. The 
function x,(H~) can be readily found from equation (4) and the table of monochromatic 
atomic continuous absorption coefficients provided by the author.!® For x,(H) the cus- 
tomary Kramers’ approximation was used as outlined, for example, in Unsdld’s book. 
Once x, = «,(H~) + «,(H) and x are both known, the ratio of the emergent flux F, to 


TABLE 2* 


ROSSELAND MEAN ABSORPTION COEFFICIENT/GRAM STELLAR MATERIAL 


| be 
10-2 1071 1 10 1? 108 
—3.97 | —2.97| -1.97| -—0.97| 0.02 | 0.96) 1.59 1.77 
| 3.97] —2.97) -1.97] -0.97] 0.02) 0.96] 1.59 | 1.89 
3.58| —2.58| —1.58| —0.61 0.57) 062 | 0.63 
—3.58| —2.58| —1:58| —0.61 0.21; 0.62) 0.95 1.62 
—3.29| —2.30| -—1.37| -—0.78| —0.61| —0.61 |......... 
—329} —2.30) -1.37) -0.78| -—0.49 0.00) 0.81 1.79 
0.8 —3.12| -106| -—0.06| 0.94 | 1.94 
—3.06| —2.25| —1.90| —2.85| 
—3.06 | -1.86| -1.58] —0.95} -0.06| 0.94 1.94 
= ~3.98 —2.96| -1.96| -0.96|} 0.04] 1.04 2.04 
—3.27| —3.02} -—2.65| -—1.91| -—0.96| 0.04 1.04 | 2.04 
10. | —4.89| —3.80| —2.89} -1.89} -0.89/ 1.41 | 2.11 
| O11) 1.11 | 2.11 
14 | —4.84| —3.84| —2.84| -1.84] -0.84| 0.16] 1.16 | 2.16 
| | 
| —4.84| —3.84| -2.84) -1.84| -0.84| 016] 1.16 | 2.16 
| —4.81| —3.81) —2.81| —1.81) —0.81| 0.19] 1.19 2.19 
} 


* Each square contains: first line, logio x (#1~); second line, logio x (H); third line, logio x. 


the flux F‘ (which is the black-body flux corresponding to the effective temperature 
7.) can at once be obtained by interpolating in Burkhardt’s table. Tables 3, 4, and 5 and 
Figures 1, 2, and 3 summarize the results of such computations. 

In Table 3 the quantities x,(4/~) per gram of neutral hydrogen and per bar of electron 
pressure and x,(/7) per gram of neutral hydrogen are tabulated for different values of the 
reciprocal effective temperature 0, = 5040/7. and for a number of selected frequencies. 


19 Op. cit., Table 2. 
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In Table 4 the final monochromatic absorption coefficients per gram of stellar mate- 
rial are given for twelve different model atmospheres; Table 4 also contains the quanti- 
ties logio F,/F{° for these models at the selected frequencies. At the Paschen and the 


TABLE 3* 
MONOCHROMATIC ABSORPTION COEFFICIENTS 


8 
vx 10-15 | ~ 
| 0.5 0.6 0.7 0.8 0.9 
1.62 2.73 4.48 714 | 109 4 
100.6 4 646 0.2207 1.070% 10-2 5.279 10-8 
2.15 3.85 6.52 10 60 1666 
55.70. | 2.573) 0.1222 5.927 10-3 (2.923% 10- 
171.6 f | 19.752 f | 10.5665 | 3.340X10-2) | 1.903K 10-3 | 
2.44 4.40 51 12.26 19 32 
130 8 7.430 0 4315 2 5441072 1518 10-3 
2.57 4 64 13.00 20.48 | 
113.2 6.432 0. 3736 2 203 10-8 1314x1073 | 
3.04 5.60 9 66 15.86 25.07 
53.26 3.026 0.1758 1.036 10-2 6.185 107! 
3.12 5.76 9.93 16.32 25.82 
38 74 2 201 0 1278 7 540 10-3 4490 107! | 
3.10 5.72 9 88 | 16 25 25 71 
29 18 1658 | 9.628x10-2| 3.388 10-4 
2.89 535 | 9.28 15 26 24 16 
17 68 1 004 3. 440 10-3 2.052107! 
2.84 5.27 9 12 15.00 23 74 
| 16.34 0.9287 0.05394 3.181 10-3 1 898% 10-4 
oe 2.79 5.17 8.95 14.73 23.31 
| 15.07 / 0.8562 0.04973 2.932 10-3 1.750% 
317.1 29 74 2811 | 0.2667 | 
en 2.55 4.74 8 22 13.55 21 48 
242.1 22.71 2.146 0. 2037 1.939 10-2 
1.00.. 2.30 4.27 7.40 12.18 19 29 
176.1 16.52 1 561 0 1481 1 410 1072 
1.36 2.53 4.39 7.23 11 46 
64.18 6 020 0 5689 5.399 10-2 5 1391073 


* Each square contains: first line, x» (H~) per gram of neutral H and for an electron pressure of 103 bars; second line, 
xy(H)/gm neutral H. The values in braces are those at either side of a discontinuity. 


Balmer discontinuities the relevant quantities are given for both sides of the discon- 
tinuity. 

Figure 1 shows a plot of the absorption coefficient for the model atmosphere for 
0. = 0.5, pe = 10% bar. 
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In Figure 2 we have plotted the values of logio F,/F\°) as a function of frequency for 
the various models. For practical purposes it is convenient to have values of logio F,/F\ 
because it is directly comparable with the observations for the gradients made on the 
continuous spectra of stars and also because it is closely related to color temperatures. 
On such a plot as Figure 2, energy distributions corresponding to any color temperature 
are approximately straight lines. Thus index marks made in this figure show the slopes 
corresponding to various reciprocal color temperatures. 

Figure 3 gives the distribution of the emergent flux in absolute units as a function 
of frequency for two of the model atmospheres considered: 0, = 0.5, p. = 10°; and 


8. = 0.8, = 10°. 


Fic. 1.—The monochromatic absorption coefficient «x» for a model stellar atmosphere with 0, = 
5040/T. = 0.5, p-= 10° bar. The ordinates denote xy in units of cm? per gram, while the abscissae denote 
frequencies in cycles per second, one small division representing 0.1 10" cycles. The ordinate of the 
horizontal line corresponds to the value of x. 


Table 5 gives computed color temperatures for the models at various frequency ranges 
of astrophysical interest, as well as the logarithm of the ratio of the fluxes on either side 
of the Balmer and the Paschen discontinuities. 


IV. COMPARISON OF THE RESULTS FOR MODEL STELLAR ATMOSPHERES 
WITH THE DATA OF OBSERVATIONS 


An important feature of the emergent flux of the radiation of stellar atmospheres is the 
sudden discontinuity at the head of the Balmer series. Since the model atmospheres show 
large variations in the amount of this discontinuity with the electron pressure, it is nat- 
ural to suppose that for a given effective temperature the average electron pressure in a 
stellar atmosphere will be indicated by the observed magnitude of the Balmer discon- 
tinuity. Now Barbier and Chalonge’’ give as the final averages of their measures of the 
quantity D = A logiy F, at the Balmer discontinuity, for more than 200 main-sequence 
stars, the following values: AO, 0.47; A2, 0.44; A3, 0.42; A5, 0.39; FO, 0.28; F2, 0.22; 
F5, 0.17; F8, 0.11. Using Kuiper’s temperature scale,’ we can refer these values for the 
different spectral types to appropriate effective temperatures; and in this manner we can 
compare the observed values of D with those given in Table 5, calculated for the model 
atmospheres. It is found that with p, = 10% we obtain values of D agreeing with those 
observed in the earlier spectral types, while a slight decrease of p, is indicated for the 
later ones. 

Barbier and Chalonge” have also measured stellar color temperatures on either side 
of the Balmer discontinuity. The spectral ranges of their measures correspond closely 
with those of the computed color temperatures B and C of Table 5. Figure 4 contains a 


20C.R., 210, 99, 1940. 
214 p.J., 88, 429, 1938. 22 Op. cit. 
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Fic. 2.—The deviations of the emergent flux Fy of model stellar atmospheres from the radiation 
F‘°) froma black surface at the effective temperature. The ordinates represent logio F,/F®) (one small 
division corresponding to 0.1 in the logarithm), while the abscissae represent frequencies (one small divi- 
sion corresponding to a frequency difference of 0.2 X 10" cycles per second). The numbers labeling the 
curves correspond to the electron pressures of the model atmospheres which they represent. The curves 
marked “H~” correspond to the results for atmospheres of pure H~. Index marks give the approximate 
slopes corresponding to the various values of @¢. 
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plot of the reciprocal color temperatures as a function of the reciprocal effective tempera- 
tures: (1) of Barbier and Chalonge’s observations; (2) of Wildt’s®* earlier calculations 
for p- = 10*; and (3) for the model atmospheres computed in this paper, also for p. = 
10*. In both cases the present models give a slope nearer that of the observed stars. It 


21/0" 


Fic. 3.—The frequency dependence of the emergent flux for two model stellar atmospheres (0, = 
0.5, pe = 108; 0, = 0.8, p = 10%). The quantity Fy (in cgs units) is plotted along the ordinate, 
while one small division of the abscissa corresponds to a frequency difference of 0.2 * 10" cycles per 
second. An auxiliary scale represents the wave length in thousands of angstroms. 


should be noted that the point with the largest value of 0, is for the sun, measured by 
Arnulf, Chalonge, and Jardin,” and is doubtless not directly comparable with the other 
observations. This is because of the fact that they may be based on slightly different 
zero points and also because the dispersions were different, thus introducing a variable 
influence of absorption lines on the gradients. Since the zero point of the Barbier and 


23.4 pJ., 93, 47, 1941. 4C.R., 210, 325, 1940. 
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Chalonge color observations may be uncertain, the vertical placing of the observed points 
in Figure 4 may be subject to uncertainties. However, it appears that the differences be- 
tween the reciprocal color temperatures for the two regions should be less subject to 
errors of the zero point than are the quantities themselves. Accordingly, in Figure 5 we 
have compared observed values of the differences in the color temperatures at Ad 4600- 
3700 and AA 3700-3150 with the results for the model atmospheres. It is seen that the 


TABLE 5* 
THEORETICAL RECIPROCAL COLOR TEMPERATURES 


O- 0.5 0.6 | 0.7 0.8 0.9 

D.. 264 0.030 

E.. 0.034 0.003 
pe=10? A.... 0.42 0.70 0.89 1.02 

59 63 72 0.84 

75 329 047 0.005 

I Seavert 0.144 0.062 0.007 0.000 
pe=108 A..... 0.36 0.54 | | 

31:: 50 68 

D 58 371 070 

(A.. 0.57 0.67 — 079 | 0.90 1.02 
Limit <B... | 51 61 68 80 0.93 

\C.. 0.48 0.54 | 0.63 | 0.73 0.83 

107 10° 108 10? 


The bottom row gives limiting values of the color temperatures for stellar atmospheres, where only 
the opacity due to H7 is effective. The electron pressure, p,*, at which this occurs, is also shown for each 
temperature. 

* Oc = 5040/Te A; B, 


10*/A = 1.40-1.80, 2.20-2.60, 2.74~-3.00, 
D = A logio F, at Balmer limit, E = A logio F, at Paschen limit. 


computed relations are not in exact agreement with the observed ones, but it is also seen 
that our present calculations represent a slight improvement over Wildt’s calculation.” 

Now the Greenwich observers have determined the color temperatures of a total of 
more than 300 stars,” in a wave-length region between 4100 A and 6500 A. This is a 
range which is nearly equal to the extremes of the ranges of the computed reciprocal color 


*® The fact that the point for the sun is in such marked disagreement with the run of the other observed 
values seems to indicate that the solar observation is not directly comparable with the others. 

26 Sir Frank Dyson, Observations of Color-Temperatures of Stars, 1926-1932, London, 1932; also M.N., 
100, 189, 1940. 
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temperatures A and B (see Table 5). Thus an approximate value of the reciprocal color 
temperature of the model atmospheres throughout the whole region may be taken to be 
the average of A and B. Figure 6 shows the relation between the reciprocal effective 
temperatures and the reciprocal color temperatures for this spectral region for (1) aver- 
ages of stars of various spectral types as observed at Greenwich; (2) Wildt’s model at- 
mospheres; and (3) the model atmospheres computed here (both [2] and |3] correspond- 
ing to p. = 10%). A value of ¢9 = 1.00 was taken for the Greenwich gradients; and so 
the whole set of observed points in Figure 6 may be subject to a slight displacement in 
the vertical direction. However, the general disagreement of slopes could not be re- 


Fic. 4.—(0, — ®,) relations as observed and computed, for the intervals \A 4600-3700 (the diagram 
to the left) and Ad 3700-3150 (the diagram to the right). The solid dots represent observations of 
Barbier and Chalonge. The single circles represent theoretical values for model atmospheres, with p, = 
103; and the double circles correspond to the earlier calculations of Wildt (Ap./., 93, 47, 1941). The 
observed point at @, = 0.72 is quite uncertain (cf. Barbier and Chalonge, Ann. d’up., Vol. 3, No. 2, 
1940). The point at 6, = 0.89 is for the sun; and, since it was obtained in a different manner from the 


others, is not strictly comparable with them. 


moved by any such translation. It is obvious that the agreement of the model atmos- 
pheres with the visual (Greenwich) gradients is less satisfactory than with either the 
photographic or the ultraviolet gradients (Barbier and Chalonge). However, the dis- 
agreements, such as they are, are all in the same sense: the 0, — 0, relation is steeper 
for the model atmospheres than that which the observations indicate. Altogether, the 
more accurate H~ absorption coefficient which has been used in the present calculations 
seems to have decreased—definitely, if not sufficiently—the disagreement between the 
model atmospheres and the true stellar atmospheres. Wildt has pointed out*’ already 
that this discrepancy cannot be traced to the effect of absorption lines, since it operates 
in the opposite direction. 

The model atmospheres show a discontinuity in F/, at the head of the Paschen series 
similar to that at the head of the Balmer series, although considerably smaller. Hall and 
Williams** have measured the quantity Dp = A log, )F, at the Paschen discontinuity for 


27 Ap.J., 93, 47, 1941. 28 A p.J., 95, 225, 1942. 
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the AO star Vega and find it to be approximately 0.092. Reference to Table 5 shows that 
for a model atmosphere with 0, = 0.5 and p, = 10% the quantity Dp = 0.104, which is 
in satisfactory agreement with the estimate of Hall and Williams. However, it should be 
remarked that Hall and Williams give their value as a minimum and state that further 
observations will be necessary to establish the exact magnitude of the discontinuity. 
Accurate investigations of the linearity of stellar gradients in the visual and the near 
infrared regions have been made recently by Williams*® and by Hall,*° respectively. 
Williams finds only slight systematic deviations from linearity in the region A 4000- 
6700, for stars of types A and F. The small differences in gradient which he does find 


a6 Q7 O08 


Fic. 5.—AQ, = 0.(AX 3700-3150) — @-(AX 4600-3700) as observed and as calculated. The ordinates 
denote A®,, and the abscissae denote 0,. The other symbols have the same meaning as in Fig. 4. 


make the color temperatures higher in the red than in the blue. Reference to Figure 2 
shows that this curvature is opposite in sign to that exhibited by the model atmospheres. 
Hall’s results for the near infrared (AX 4600-10,000) support the assumption that the 
main-sequence stars of types A and F radiate essentially as black bodies in this region of 
the spectrum, except for the Paschen discontinuity. On the scale of Figure 2 the curva- 
ture found by either Hall or Williams would not be noticeable over the ranges of their 
respective observations. 

To summarize: The model atmospheres studied in this paper, with the help of the 
continuous absorption coefficient of H~ as recently recomputed, predict an effective tem- 
perature to color temperature relation which fits observed curves fairly satisfactorily in 
the ultraviolet, less so in the photographic region, and still less so in the visual region. 


29 Pub. Obs. U. Michigan, 7, 93, 147, 159, 1939. 30 Ap.J., 98, 231, 1942. 
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In all three cases, however, the theoretical relations derived here lie somewhat closer to 
the observational ones than those obtained earlier by Wildt, who used a less accurate 
absorption coefficient for H-. 

The models treated in the present paper show a definite trend of the color temperature 
to increase with increasing frequency, while actual stars seem to show a fairly constant 
color temperature, which, if anything, decreases slightly with frequency. 
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Fic. 6.—The (6,-9,) relation in the spectral region \A 6500-4100, as observed and as calculated. The 
filled circles represent the Greenwich observations. The single circles correspond to the values calculated 
in the present paper with p, = 10%, while the double circles represent the values according to Wildt’s 
earlier calculation. 
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friendly encouragement, without which this paper would have been impossible. Dr. Jesse 
L. Greenstein is to be thanked for helpful discussions of observational data on stellar 
spectrophotometry. The writer’s wife, Dorothy C. Williamson, was of great assistance 
in much of the computing and in typing the manuscript. 
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OBSERVATIONS OF THE LIGHT OF THE NIGHT SKY 
WITH A PHOTOELECTRIC PHOTOMETER* 


C. T. ELVEY 


ABSTRACT 


A preliminary report is given of the observations of the brightness of the night sky with a photoelectric 
photometer sensitive to the infrared region of the spectrum. A brief description of the apparatus is 
given. Curves of the nocturnal variation of intensity of the light of the night sky have been obtained on 
45 nights, and it is found that they may be classed in two groups: one with a minimum near midnight 
and the other in which the intensity decreases through the night. 

The irregular variations of intensity during the night, which are greater and more numerous when 
there is considerable magnetic activity, are not the same in different directions in the sky, thus showing 
the nonuniform character of the light of the night sky at such times. 

Tests to determine the height in the atmosphere and the motions of the irregular patches of light 
have given no positive results. Also no correlations were found between the fluctuations in the light of 
the night sky and the variations of the ionosphere as indicated by the signal strength of reflected radio 
waves, the night sky being observed around the downcoming radio wave. 


Many interesting results have been obtained from spectrophotometric observations 
of the light of the night sky; and, although this method has certain distinct advantages, 
there are also some disadvantages: the accuracy is limited and the length of time that 
it takes to obtain a single observation is long. The latter is a particularly important 
disadvantage in connection with studies of the rapid changes in brightness which occur 
in the light of the night sky. The photoelectric photometer is valuable for this work, 
since it combines accuracy with a rapid indication of the intensity, a single observation 
taking only a few seconds or perhaps two or three minutes when the greatest sensitivity 
of the instrument is required. 

This paper is intended as a preliminary report on the photometry of the light of the 
night sky with a photoelectric photometer. It is a brief presentation of the observations 
which have been made during the past year. Several problems have arisen, but their 
solutions will have to wait until much more observational material can be collected. For 
instance, additional observations will be needed to redetermine the variation with zenith 
distance of intensity of the light of the night sky in the infrared region, so that the height 
of the auroral layer may be checked. The preliminary observations with the photoelec- 
tric photometer indicated a height of 125 km for this layer,! whereas the spectrophoto- 
metric observations of Elvey and Farnsworth? for the stronger lines in the red-green 
regions of the spectrum indicate a height of 500 km. 

Three photometers have been constructed for this work, all of them having the same 
general design. Slight variations, however, were made both in the optical systems and 
in the amplifiers. In one photometer the entire unit excepting a storage battery is 
housed in a box. In the others the optical system, the photoelectric cell, and the ampli- 
fier tubes were separated from the control cabinet by a shielded cable, thus making it 
possible to operate the instrument from a lighted room. 

The optical system consists of a condensing lens whose aperture is 109.5 mm and 
whose focal length is 139.7 mm. The lens focuses the sky on an aperture of a metallic 
box which shields the electrical system. The photoelectric cell is placed immediately 


* Contributions from the McDonald Observatory, University of Texas, No. 62. 


' A paper presented at the dedication of the National Astrophysical Observatory at Tonanzintla, 
Puebla, Mexico. 
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behind the aperture; and, in order that it may be effectively shielded, a copper screen 
is soldered over the aperture. This reduces the intensity of the light slightly, but it 
greatly increases the stability of the electrical circuit. A sharp-cutting red filter, Corning 
No. 241, removes most of the visual part of the spectrum and only that part of the spec- 
trum which is greater than \ 6400 is effective. A shutter operated by a bulb and a long 
rubber tube serve to expose the photoelectric cell to the night sky. 

The receiver is a caesium oxide photoelectric cell, No. 71-TA, made by the G-M Lab- 
oratories. Its maximum sensitivity is \ 7000. At ordinary temperatures these photo- 
electric cells have a dark current—mostly thermal emission from the cathode—which is 
five to ten times the size of the photoelectric current when the cell is exposed to the 
night sky. Slight temperature changes will thus vitiate the observations. A gust of 
wind causes fluctuations of the indicating meter which are sometimes larger than that 
produced by the photoelectric current. It was found by experiment that the tempera- 
ture variations could be sufficiently dampened by surrounding the photoelectric cell 
with a laver of cork an inch or more in thickness and by using a double window to admit 
light to the cell. This produces a steady drift in the electrical circuit for several hours 
at the beginning of the night when the temperature has dropped several degrees, but it 
is slow enough that the zero may be adjusted and read only once prior to each observa- 
tion. The observations are made at 5-minute intervals throughout the night; but during 
twilight and dawn, when the intensity is changing rapidly, the rate of making readings 
is doubled. 

The photoelectric current is measured by means of a direct-current amplifier of the 
negative feedback type described by Roberts.* However, some slight changes and ad- 
justments were made, since the amplifier is always used at its highest gain. The current 
is read on a microammeter with a scale of 0-50. The power supply is from batteries. 
Two of the instruments use storage batteries for the filament current, and the other 
uses low-drain tubes and employs dry cells for the filament current. The other poten- 
tials are from small B batteries. 

Each photometer is equipped with a standard lamp built into the lid which covers 
the lens, thus making it possible to test the instrument in a lighted room. The intensity 
of the light is adjusted so that it produces a deflection of the same order or magnitude 
as that from the night sky. The lamps are similar and can be readily removed so that 
intercomparisons can be made. One lamp has been set aside and designated as a pri- 
mary standard and all the other lamps are compared with it. The standard lamp has 
given a check on the stability of the instruments through the night and from night to 
night. If the problem can be continued, it is expected that the lamps will give a means 
of intercomparing the brightness of the night sky as observed on different nights and 
at different localities. 

It is possible to observe the sky with some moonlight, since this photometer is sensi- 
tive to the infrared portion of the spectrum. Tests have been made on nights of photo- 
metric quality with the full moon in the sky. The brightness of the region around 
the pole was increased by only 1 mag. over a night without a moon. Whenever the sky 
is free from haze and thin clouds, it is possible to observe until the moon reaches quarter- 
phase, with hardly any noticeable influence on the brightness of the sky. With moon- 
light there is no trouble in distinguishing whether haze is present, for then the sky is 
abnormally bright. Also a very thin cloud passing across the field of the photometer will 
cause a large, sudden change in brightness. On the other hand, the haze and thin clouds 
without a moon in the sky have very little effect upon the brightness of the sky. There 
is merely a smoothing of the irregular variations, since the haze integrates the light from 
the entire sky. 

During the 14 months ending in June, 1942, the curves of the nocturnal variation of 
intensity of the night sky have been obtained on 45 nights with only a few of them 


3 Rev. Sci. Inst., 10, 181, 1939. 
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having 1 or 2 hours missing. In addition, many parts of nights have been observed. 
For all these observations the instrument was pointed toward the north pole. 

An inspection of the general shape of the curves indicates that they may be grouped 
into two classes, one which has a tendency to produce a minimum of light around local 
midnight or shortly after and the other which shows a tendency for a steady decrease 
in the light through the night. In both cases there are irregular variations, but in the 
latter they seem to be more pronounced. There are a few cases when it is rather uncer- 
tain to which class the curve might be assigned. It is further noticed that there is a 
tendency for a grouping of the curves of a given class. For two or three nights the curves 
would be similar. During this period of observation there was found no evidence that 
the classes of curves are related to the seasons; however, the examples cited below hap- 
pen to have one group all of which are winter observations and the other all spring or 
summer observations. 

Examples of the class of curves for which there is a minimum of intensity in the 
vicinity of local midnight are shown in Figure 1. The dates and abscissae are expressed 
in Universal Time. Some of the curves are quite free of the irregular fluctuations in in- 
tensity, while others have moderately large ones. The last curve in the diagram is one 
of the few which could not be assigned to either of the two classes. 

In Figure 2 are presented a series of curves in which there is a general decrease in the 
intensity of the light of the night sky through the night. Different degrees of the irregu- 
lar fluctuations are shown. It has been known for some time that the irregular fluctua- 
tions are associated with times when the terrestrial magnetic activity is great; and re- 
cently D. N. Barber‘ has shown an excellent correlation between the irregularities of 
the light of the night sky as observed at the Lick Observatory and the terrestrial mag- 
netic activity as observed at the Mount Wilson Observatory. 

The observations have been extended into the evening and the morning twilights; 
and, since the region observed is that of the pole, it is symmetrically placed with respect 
to the setting and the rising sun. Since the photometer integrates both the continuous 
and the emission spectra over a large range of wave lengths, we cannot be certain 
whether the twilight effect is a result of scattered sunlight or an emission spectrum. By 
comparing the morning and evening observations at the pole we find that the change 
of intensity is symmetrical with respect to the position of the sun. If it were caused by 
an emission spectrum, we should expect differences between the morning and the evening 
intensities. For the evening twilight the sky will have been in the sunlight for the pre- 
ceding hours, while for the morning twilight it will have been in the shadow for a similar 
length of time. We can conclude that the observed enhancement at twilight is a result 
of scattered light from the twilight portions of the atmosphere. 

Observations have been made with two photometers pointed to different parts of the 
sky, and these have shown that the irregular fluctuations of intensity are not the same. 
Hence it must be concluded that the light of the night sky is not uniform as observed 
from one point and that the irregularities are probably due to patches of light. Since it 
was thought that the light from the night sky comes from the region of the £-layer,! 
a series of observations have been planned to determine the character of the irregular 
fluctuations. These have not been completed, but the observations of the following 
paragraphs indicate that the luminous areas are probably not in the E-layer; the series 
of observations will have to be greatly extended. 

Observations have been made with one photometer at the McDonald Observatory 
pointing toward the pole, while another instrument was taken to various positions in 
the near vicinity. The purpose was to determine the height and the motions of the 
patches of light which produced the irregular variations during the night. One of the 
first sets of observations indicated a lag of several minutes between the fluctuations of 
intensity as observed at Mount Locke and at a point near Fort Davis, a distance of 


‘ Lick Obs. Bull., 19, 105, 1941. 
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12 miles approximately in the east-west direction. However, repeated trials have not 
verified this result. The longest distance that has separated the two stations to date is 
approximately 25 miles. If the problem is to be carried on, the stations should be set 
up at least 100 miles apart. 
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Fic. 1.—Curves of nocturnal variation in the intensity of the light of the night sky, showing minima. 
The dates of the curves are February 18 and 20 and March 17, 1942; November 18 and December 14, 


1941. 


An attempt was made to obtain observations over a base line of 425 miles, between 
McDonald Observatory and Forth Worth, Texas, but poor observing conditions have 
made it impossible to obtain simultaneous observations. The observations at Fort 
Worth were being made through the kind co-operation of Mr. Dan C. Taulman, an 
amateur astronomer, but have had to be discontinued because his work in defense proj- 
ects now consumes most of his time. 
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Fic. 2.—Curves of nocturnal variation in the intensity of the light of the night sky, showing a decrease 
of intensity through the night. The dates of the curves are as follows: May 9, 17, and 18, 1942; June 27 
and July 22, 1941. 
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The agreement between observations taken with two photometers has been checked 
a number of times. The bottom diagram of Figure 2 shows the agreement between ob- 
servations taken at McDonald Observatory (Mount Locke) and a point near Fort 
Davis, a distance of 12 miles approximately east-southeast. These observations were 
made on a night in which there was some large variation in the intensity of the light 
of the night sky; it will be noticed that most of the features are quite accurately repro- 
duced by the two instruments. There are a few discrepancies which are well above the 
accuracy of the instruments, showing that for short durations the light was of different 
intensities as observed from the two locations. Notice particularly the difference at 
5h 20". 

On this night I was operating the instrument at Fort Davis and had excellent oppor- 
tunities of watching the sky. The readings of the meter were made with a faint flashlight 
and my eyes were dark adapted. The brightness of the sky changed greatly in the 
infrared region, but the eye had difiiculty in noticing that the sky was brighter at one 
time than at another. In my notes I have recorded that at 7" 00™ the northern horizon 
appeared brighter than the southern horizon. 

The irregularities in the light of the night sky probably represent areas of increased 
brightness over that of the background, and it is reasonable to suppose that these areas 
are regions of greater excitation. If the height to the layers in which the radiation is 
produced is that indicated by the preliminary observations, i.e., the £-layer, then we 
might expect to find correlations between the irregular fluctuations in the light of the 
night sky and the characteristics of the ionosphere. Bradbury and, Sumerlin,’ observing 
at the Lick Observatory with a photoelectric photometer and making observations of 
the ionosphere at Stanford University, in California, did not find any correlation between 
the brightness of the sky and the ionospheric conditions; however, it is possible that the 
disturbed areas of the upper atmosphere are too localized to give a correlation between 
these two localities. We do not have sufiicient equipment available to conduct investi- 
gations of the ionosphere with pulse transmitters simultaneously with our observations 
of the brightness of the night sky. However, it might be expected that the radio signals 
reflected from the ionosphere would give some idea of the variations in the ionized 
layers; and we have tried observing simultaneously the signal strength of a radio station 
and the brightness of the night sky in the region surrounding the computed position of 
the downcoming radio wave. The field of view of the photometer is 8°, and it is very 
probable that this is sufficiently large to include the difference between the computed 
and the actual positions of the downcoming wave for signals in the broadcast band, but 
this may not be the case for the shorter waves which are reflected from the /’-layer. For 
the latter case the position of the photometer has been shifted from night to night. 

The preliminary trials were conducted by using the broadcast station KRLD, in 
Fort Worth, Texas, whose frequency is 1080 kilocycles. There was a suggestion of a 
correlation of the frequency of the fading of the signal with the intensity of the light of 
the night sky, but there were no correlations with signal strength. 

Further tests were tried by using the signals from radio station KRCA in San Fran- 
cisco, California, which is broadcasting on a frequency of 9.48 megacycles. Some excel- 
lent records were obtained of both the radio-signal strength and the brightness of the 
night sky; but no correlations were evident, although both showed well-defined varia- 
tions. 

The negative results given above and also those of Bradbury and Sumerlin® point 
toward the origin of the radiation as being in the high atmosphere, in contradiction to 
the preliminary observations of height as given by the photoelectric photometer, per- 
haps from the same region as the green auroral line, about 500 km above the surface of 
the earth. 


5 Terr. Mag. and Elec., 45, 19, 1940. 
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The time constant of the electrical circuit of the photometer is an appreciable fraction 
of a minute, and consequently variations of intensity of the light of the night sky which 
might be of the order of a second or so would hardly be detectable. During the times 
when the instrument was exposed to the sky at or near the maxima of the irregular varia- 
tions, it was noticed frequently that the needle of the indicating meter was unsteady. 
Tests have shown that this unsteady condition is not of instrumental origin. With no 
light on the cell the needle of the meter is steady, and this is also the case if the cell is 
exposed to the standard lamp. This would indicate that there are very rapid fluctua- 
tions of intensity of the light of the night sky, but this should be verified with a photom- 
eter which is constructed to detect rapid changes of intensity. 


This work has been greatly assisted through grants of financial aid given both by the 
American Academy of Arts and Science and by the American Philosophical Society, for 
which I am very grateful. It is a pleasure to acknowledge also the assistance of Mrs. C. T. 
Elvey and of Mr. Joseph Ross Brown, each of whom helped to obtain long series of ob- 


servations. 


McDOoNnaLp OBSERVATORY 
September 1942 


TENTATIVE IDENTIFICATION OF THE HERZBERG BANDS OF 0, 
IN THE ULTRAVIOLET SPECTRUM OF THE NIGHT SKY* 


P. SwIncs 


As is well known, the identifications of the night-sky bands are still very unsatisfac- 
tory; even the strongest ultraviolet feature at \ 3556, which appears always in long 
spectroscopic exposures on faint nebulae, is still unidentified. Some time ago, I had the 
privilege of receiving from Dr. J. Dufay, of Lyons (France), a manuscript copy of an 
interesting paper in which a rather convincing identification is suggested for the strong 
ultraviolet bands, including \ 3556.! Dufay believes that these may be due to the “‘for- 
bidden”’ transition — (Herzberg system) of O2. This band system has been 
observed only in absorption in the region \ 2400-A 2900. The absorption consists of a 
progression v’’ = 0, each (v’, 0) band having only a Q branch (AK = 0) degraded to 
the red; the lower electronic level is the ground state of Oz. Dufay suggests that several 
night-sky emission bands are the (v’ — v’’ > 0) transitions of the Herzberg system. 

The v” = 0 progression, which was found in laboratory absorption by G. Herzberg,’ 
has been measured by several authors’ either in the laboratory or in atmospheric absorp- 
tion. The measured wave lengths of the origins agree within about one angstrom unit 
for the bands of \ < 2637. Herman has measured three additional weak bands to the 
red of \ 2637, at AX 2685.0, 2737.6, and 2795.2; and there may still be two more at 
approximately \ 2843 and \ 2912. Hence the actual v’ designations are still in doubt, 
but the energies F’ (v’) are fairly well known. 

To obtain the v’ = constant progressions, Dufay determines the v(v’, 0’’) formula of 
the system by least-squares solutions, accepting for the v”’ term a simplified expression 
obtained from the Schumann-Runge and the atmospheric O2 bands, which have the 
same lower electronic level *2>. Then Dufay compares the calculated values of X(2’, v’’) 
with his wave lengths of the night-sky features and finds several cases in satisfactory 
agreement. 

I have undertaken a similar discussion on the basis of the results on the night-sky 
spectrum obtained at the McDonald Observatory.‘ Instead of using least-squares solu- 
tions, I have simply taken for the upper state the “best” »(v’) corresponding to the 
observed Herzberg absorption bands.° For the ground electronic state I have used the 
formula, including the v’’* and v’’4 terms, obtained by J. Curry and G. Herzberg,® which 
is better for large values of v’ than the simplified quadratic formula used by Dufay. 
More refined calculations are useless. The wave lengths of the night-sky spectrum are 
not accurately known, and we do not have accurate information on the widths or pro- 
files of the bands. Moreover, we do not know what the rotational-intensity distribution 
should be in the night-sky bands of a specific molecule. From the observations we 
gather that ihe rotational-intensity distribution is of a low-temperature type. Yet we 


* Contributions of the McDonald Observatory, University of Texas, No. 63. 

! Dufay’s note has probably been published in the Comptes rendus, Paris. See Bull. 3 of the Committee 
for the Continued Distribution of Astronomical Literature, October, 1941. 

2 Naturwissenschaften, 20, 577, 1932. 

3 Chalonge and Vassy, C.R., 198, 1318, 1934; Herman, Ann. de Phys. (11th ser.), 11, 548, 1939; 
Miss N. Morguleff and Mrs. A. Vassy, Ann. d. Ap., 1, 427, 1939. 

4 Elvey, Swings, and Linke, Ap. J., 93, 337, 1941. 

5 The uncertain bands \ 2843 and \ 2912 are not retained; hence v’ = 0 is adopted for \ 2795.2. 

6 Ann. der Phys., 19, 800, 1934. 
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cannot, a priori, determine the amount by which the observed rotational centers of 
gravity of specific bands should be shifted to the red, relative to the band origins. Such 
a temperature shift may even be different for different values of v'!’ 

Taking into account the various sources of uncertainty, a number of coincidences are 


found, the majority of which are probably real. They are listed in Table 1: column 1 


TABLE 1 
HERZBERG BANDS OF O2 IN THE NIGHT-SKY SPECTRUM 


dX AND INT. IN 


| 

| 

dX AND INT. IN | 
NIGHT SKY | 
| 


Oz IDENTIFICATION i O2 IDENTIFICATION 
| 
| 


NIGHT SKY 


} 
Transition d (Calc.) | Transition | A (Calc.) 
(6,5) | 3111 |} 3598@).......... ((6,8) | 3585) 
} 
3141 || 3636(2).......... (1,6) | 3631 
| (3,4) 3143 (3,7) | 3637 
3157 3164) || 3664 (1)... (5,8) | 3656 
| 
3205 || 3707(3).......... | (7,9) | 3702 
(0,3) 3208 | | 
(2,4) 3209 || 3742 (4).......... (0,6) | 3733) 
| (4,8) 3737 
3233 (0) 3224 | 
3263 (1).......--. (6,6) 3257 | (38) | 3831 
| 
3298 (3)... (3,5) 3294 3948 (3)......... (4,9) 3939 
| (0,7) 3942 
| 3984(2)..........| (6,10) | 3973) 
(2,5) 3367 || | 
(0,4) 3368 4048 (2)..........) (1,8) 4044 
(4,6) 3381 | (3,9) | 4044 
3051). 3414) 4072 (5). (5,10) | 4060) 
3460 (2) (1,5) 3451 || 4174 (5)..........) (4,10) 4160) 
(3.6) 3459 || (0,8) 4172 
| 
3488 (3) } (5,7) 3478 | 
3556 (10nn) (20) 3539) | 
| (05) 3542 | | 
| (4,7) 3552 | 
| | 


contains the night-sky wave lengths and intensities, as found at the McDonald Observa- 
tory;! columns 2 and 3 contain the identifications which should be added to the ones 
previously given.‘ The v’ values are uncertain, and they may possibly have to be in- 
creased by one or two units. The wave lengths in column 3 refer to the band origins, 
which should lie slightly to the violet of the observed maxima of the night-sky lines. 
The O: bands provide probably the major contributions for all the lines of wave length 
\ <3598, excepting possibly \ 3378 and \ 3425 but including the strong, unidentified 
lines. For \ > 3598 the influence of the Vegard-Kaplan system becomes important. 


7 This shows again the importance of renewed attempts to obtain the ultraviolet spectrum of the night 
sky with a better resolving power than has hitherto been available. 
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Assuming that the identification is correct, the average red shift of the observed . 
bands, relative to the calculated band origins, corresponds to a temperature of the order 
of 150° K. 

In a number of details the identifications given in Table 1 differ from the ones sug- 
gested by Dufay, but they favor his general hypothesis. As mentioned by Dufay, the 
excitation potential of the Herzberg bands is of the order of 4.5 volts; hence, it is not 
very different from the energy of the upper level of the green [O 1] line (4.2 volts). 
The emission of the O2 bands may be expected if the luminescence is excited according 
to Chapman’s mechanism (recombination of two oxygen atoms). According to Herz- 
berg,® the °={ — 827 system is more likely to appear than the Schumann-Runge bands 
of Or. 


McDONALD OBSERVATORY 
August 1942 


8 Private communication. 
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RECENT PROGRESS IN ASTROPHYSICS' 


OBSERVATIONS OF THE PLANETS BY LYOT, GENTILI, AND 
CAMICHEL FROM THE PIC DU MIDI IN 1941 AND 1942 


Letters and photographs have been received in this country from MM. Bernard Lyot 
and Marcel Gentili describing numerous excellent observations of the planets made by 
them and M. Camichel at the Pic du Midi in 1941 and 1942. Prior to 1941 the largest 
telescope at the Pic had been a 23-cm. refractor. For the 1941 opposition of Mars a 
38-cm. objective was borrowed from the Observatory of Toulouse and installed at the 
mountain station. The observing conditions were evidently very superior, and the lens 
proved to have excellent definition. According to Lyot, “it showed finer detail on the 
planet (Mars) than had been observed with the 83-centimeter refractor at Meudon under 
the best seeing conditions.” 

The 1941 opposition of Mars was studied visually and photographically at the Pic 
du Midi. The visual observations are combined and presented in the accompanying 
map of Mars. At first the observers felt that some of the details were too geometric to 
be real, but the canals and other fine details of the map have been fully confirmed by 
photographs made with the same instrument. 

The several hundred photographs of Mars which were secured during the 1941 oppo- 
sition were examined and studied at Meudon the following winter. In particular, Lyot 
prepared a series of 18 composite photographs of the planet. For each of these positives, 
several negatives were superimposed, thereby eliminating the chance grain or dust 
effects of a single photograph but enhancing the contrast for the faint details common to 
all of them. These composite pictures are on glass and were shown at a meeting of the 
Société Astronomique de France early in May, 1942. No copies have as yet been re- 
ceived in this country. The photographs of Mars (Pls. II and III) are reproductions of 
individual pictures made at the Pic du Midi in September and October, 1941. 

In addition to these observations of Mars, studies have been made of Mercury, Sat- 
urn, Jupiter, and Jupiter’s satellites. During the summer of 1942, Mercury was studied 
photographically and visually. Numerous details were seen on the little disk, and the work 
was continued until the planet was within a degree and a half of the solar limb. Photo- 
graphs of Saturn show a faint band and a small bright spot on the surface of the planet, 
and the disk is clearly visible through the outer ring. Surface markings on the four 
bright satellites of Jupiter were carefully examined. Each observer made an independent 
series of drawings over many nights, and the mutually confirmed details proved to be 
sufficient for the preparation of maps for each of the satellites. 

HELEN W. Dopson 

WHITIN OBSERVATORY 
WELLESLEY, MASSACHUSETTS 

November 1942 


! |This new section of the Journal will be devoted to critical summaries of research carried out abroad, 
which, because of war conditions, is not generally available to American astronomers. In some cases 
proofs or single copies of important foreign publications have reached this country; in others the informa- 
tion will be based upon private communications from the authors or from astronomers in neutral coun- 
tries. In some instances the original authors have expressed the wish that their work be recorded in an 
American journal. Since in most cases it is not possible to send manuscripts or proofs to the original 
authors, the responsibility for these articles will rest with the American reviewers, and their names will 
appear at the end of each article —Ep1Tor. | 
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NOTE ON DR. STERNE’S REVIEW OF VOLUME VI OF THE 
ANNALS OF THE ASTROPHYSICAL OBSERVATORY 
OF THE SMITHSONIAN INSTITUTION 


The Editor has kindly sent me advance proof of Dr. Sterne’s review of Volume VI 
of the Annals of the Smithsonian Astrophysical Observatory. It seems unfair to my 
colleagues and to the readers of the Journal to leave them with the prevailing impression 
that our campaign of measurement of the solar constant of radiation amounts only to 
this: it is acknowledged to be a worthy example of perseverance and accurate measure- 
ment and may be compared to a large pearl—costly, beautiful, worthy to be preserved, 
to be admired, but useless. 

I therefore quote from Dr. Sterne as follows in regard to the variation of the sun: 
“The crucial test . . . . is to see whether independently obtained measures, at different 
stations show significant correlation.”’ Since receiving the advance proof of his review, 
I have made such a crucial test by comparing cases of well-marked day-to-day rise and 
fall of solar radiation, as given in the Annals and observed at Montezuma, Chile, with 
corresponding cases of rise and fall of solar activity of another kind, universally be- 
lieved in, reported from an observatory of the Old World. The cases selected came 
from occurrences in two autumn months scattered throughout the years 1923-1937. 
The comparisons were 75 in number, and the average change of the solar constant in 
them, as measured at Montezuma, was 0.7 per cent. The correlation coefficient between 
these wholly independent variables is 43.1 + 6.3 per cent. 

As one of the independent variables is undoubtedly a solar phenomenon, ordinarily 
observed on the sun’s surface by photography, and as the coefficient of correlation is 
seven times its probable error, I feel that only prejudice will deny that this comparison 
indicates the reality of day-to-day solar-constant variation, as measured by Smithsonian 
observers. 

I cannot agree with Dr. Sterne as to the absence of correlation between the variation 
of the solar constant and the variation of Greenwich faculae areas, as plotted in Figure 17 
of the Annals. For first having removed the average solar constant and the average 
facula number for each solar rotation separately from the individual values plotted in 
Figure 17, so as to leave the variations unaffected by the march of sunspot activity 
and of solar radiation periodicities, I computed the correlation coefiicient between the 
remaining departures of the two variables. From 173 pairs the result is: 


r = 39.3 + 4.3 per cent. 


The correlation coefficient is thus nine times its probable error. If I had computed the 
correlation with regard to the fact that the solar-constant changes anticipate the 
changes of facula numbers, as suggested at page 196 of the Annals, the correlation would 
have resulted still higher. 

I will not prolong this note to prove that Dr. Sterne is also in error in stating that 
the solar rotation period is not shown by our observations of the solar constant. 

As to the reality of the long-period solar changes, whose unorthodox derivation Dr. 
Sterne condemns, I may recall that he himself verified some of them from our prelimi- 
nary publication of the observations.' Using his own approved mathematical method 
of analysis, he found probabilities of the reality of seven of them, ranging from 10° to 
10'°. It may be a fatuous satisfaction; but I like to see, from Figure 14 of Volume VI 
of the Annals, that the summation of these unorthodox periodicities represents the 


1 Proc, Nat. Acad., 25, 562, 1939. 
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curve of original monthly means from 1920 to 1938 with great fidelity. I look forward 
hopefully to the year 1945 for a crucial test of their excellence, arising from our predic- 
tion through them of a large decrease then in the intensity of solar radiation. 

Again applying Dr. Sterne’s criterion, quoted above, to verify the long solar periodici- 
ties, I will refer to my paper ‘‘An Important Weather Element Hitherto Generally Dis- 
regarded.’” Here, in Tables 1, 2, and 3 and Figure 7, it is shown by 37 repetitions that 
our solar period of 8{ months has continued in the temperature records of Vienna, 
Copenhagen, and New Haven for 140 years unchanged in phase except as affected by 
yearly seasonal influences. The temperature range is 1.15° C. I have in manuscript the 
proofs that others of these solar periods have also affected the weather to a controlling 
degree during that long interval. This evidence is based upon the records maintained 
by the official weather bureaus of these several countries, which are surely sources inde- 
pendent of the Smithsonian Institution. 

Thus I feel encouraged in the belief that the Smithsonian campaign of solar-constant 
measurement is not merely a useless, costly jewel of exact measurement, to be looked 
upon as a classic gallery specimen, as critics would imply. I can see that it is destined 
to furnish to meteorology an important and, indeed, controlling weather element, hither- 
to unused. I regret that I am at a disadvantage in making this rejoinder, because at 
this time it would be improper to publish my evidence in full. Should Dr. Sterne be in 


Washington, I will be glad to show it to him confidentially. 
C. G. ABBot 


WASHINGTON, D.C. 
November 1942 


CORRECTED POSITION AND SPECTRAL CLASS OF FURUHJELM 31 


A curious error has come to light as a result of the accidental discovery of a dM2 star 
of the eleventh photographic magnitude on the McCormick spectroscopic survey plates. 
The star is AC+44°472.15 (a= 11"5™50219, 6= +43°58'24"0)! with a proper motion of 
0753 found by Furuhjelm and published as No. 681 in his ‘Recherches sur les mouve- 
ments propres ...., I’ and again as No. 31 in a reprint* in which large proper motions 
only were listed. However, both publications contain a clerical error of 2° in the declina- 
tion of the star. That the declination should be +43° and not +45° is easily ascertained 
by consulting the references in column 8 of Furuhjelm’s catalogue.” 

Unfortunately, the error has been carried to several other publications. It is repeated 
in Luyten’s “List of Stars... .”’4 (star No. 293), in the Cincinnati Catalogue’ (star No. 
615), and in Kuiper’s ‘‘Notes on Proper Motion Stars.’’® As far as could be found from 
the literature, the only attempt to observe the star was made by Kuiper. He observed 
a star of 12.5 visual magnitude at the erroneous position and determined its spectral 
class as K1, from which he concluded that it must be regarded as a subdwarf. As might 
be expected, under the circumstances no trigonometric parallax has been determined as 
yet. 

A. N. VyssoTsky 
LEANDER McCormick OBSERVATORY 
December 1942 


2 Smithsonian Misc. Coll., 101, 1941. 

1 Helsingfors Astr. Catalogue, 4, 1903. 

2 Acta Soc. Sc. Fennicae, Tome 48, No. 1, 1916. 

3 Finska Vetens.-Soc. Foérh., 59, 1916-17; Afd. AN:o, 22, 1917. 

4 Lick Obs. Bull., 11, 1, 1923. 

5 Pub. Cincinnati Obs., 20, 1930. 6 Ap. J., 92, 126, 1940. 
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Magnitudes and Colors of Stars North of +-80°. By F. H. SEAREs, F. E. Ross, and M.C. Joyner. 
(“Papers of the Mount Wilson Observatory,” Vol. VI; Carnegie Institution of Washington 
Publication 532.) Baltimore: William Byrd Press, 1941. Pp. 89+ iii. Paper bound, $1.50; 
cloth bound, $2.00. 


The Mount Wilson Catalogue of magnitudes and colors for 2271 stars surrounding the North 
Pole is the fruit of nine years’ photometric work and has taken its place as one of the few collec- 
tions of fundamental photometric standards. 

The requirements of photometric standards, though simple to state, are difficult to meet. 
Only those who have had practical experience in photographic photometry are aware of the 
severe requirements by which the three essentials (accuracy of scale, determination of zero 
point, and knowledge of color coefficient) are circumscribed and may be attained. The difficulty 
of the problem may be measured by the rarity with which these requirements are met. 

The difficulty of the problem is equaled by its importance. Few branches of modern astron- 
omy are independent of the magnitude scale, though the relationship is not always immediately 
apparent. In determinations of the distances of stars or stellar systems the apparent magnitude 
is as important a factor as the period of a Cepheid; discussions of absolute magnitudes of groups 
of stars on the basis of their motions depends no less on the apparent magnitudes: to give a 
concrete example, the apparent magnitudes of most red variable stars that have been based 
on the B.D. scale (and their great number is apparent only to the student of the subject) will 
be revised downward by about a magnitude when accurate magnitudes for these stars are de- 
termined—and their absolute magnitudes will likewise be diminished. 

Standard photometry, then, is not a subject to be studied only for its own sake. It is an 
urgent necessity for all other branches of the science. Dr. Seares and his co-workers have es- 
tablished the backbone of a reliable and consistent system, first in his work on the North Polar 
Sequence and, now, in the extension of the standards that is contained in the volume under 
review. The importance, and the ungratefulness, of the task may both be measured by the 
remark in the Introduction: ‘‘Ten years after their adoption [by the I.A.U.] the relation of the 
various photometric catalogues to the NPS system was still mostly unknown and some of the 
current investigations gave little evidence that any such system existed.” (Italics the reviewer’s.) 

At its last meeting the International Astronomical Union appointed a subcommittee to deal 
with problems of standard sequences, and it is greatly to be regretted that the labors of such a 
committee should be in abeyance. Four principal needs in the domain of standard sequences 
should be focused by such a committee: (1) the need for standard magnitudes between the 
thirteenth and twentieth magnitude, south of the equator; (2) standards of photovisual magni- 
tudes are needed for faint stars, especially south of the equator; (3) several groups of carefully 
standardized magnitudes are desirable in intermediate declinations, for use in investigations 
where the North Polar Sequence cannot be used; (4) a catalogue raisonné of existing published 
sequences. 

The Mount Wilson studies of standard magnitudes are essential bases for the first three of 
these four requirements: (1) will rely on the Mount Wilson photographic magnitudes in Selected 
Areas, (2) would be greatly indebted to the catalogue now under review, and (3) must be based en- 
tirely upon that catalogue. There has been some question in recent years concerning the desira- 
bility of keeping to the International System of Magnitudes as originally defined, partly because 
the aluminization of mirrors has changed the color system, partly because the stars at the North 
Pole itself, used to establish the zero point, are found to be obscured and therefore not of nor- 
mal color for their spectral classes. But the system has established its usefulness by the most 
important of all tests—accuracy and consistency—and these have made it the accepted system 
for pragmatic reasons. The Mount Wilson Catalogue of stars north of +80°, like the North 
Polar Sequence and the Mount Wilson Catalogue of Selected Areas, is a cornerstone of present- 
day astronomy. As was said of the work of an eminent physicist, ‘‘When other people publish 
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results, we say that they should be checked and confirmed; but when he publishes something, 
we say it is so.”’ It is doubtful whether the photographic standards surrounding the North Pole 
will be redetermined in the near future. 

CecILiA H. PAYNE POSCHKIN 


Harvard Observatory 


Essentials of Astronomy. By JoHN CHARLES DuncAN. New York: Harper & Bros., 1942. 

Pp. 181. $1.85. 

This new text on astronomy impresses the reader as being different in several ways from 
the usual college text. The author has used a novel arrangement of the material, presenting 
first ‘‘a description of the sky as it appears to the most casual observer,’’ then ‘‘the changes in 
that appearance that take place with changes in the observer’s position and with the lapse of 
time,” and the interpretation of “these appearances as to the simple and majestic motions of 
the heavenly bodies.’”’ The author goes on to present Kepler’s laws, the law of gravitation, and 
the physical laws underlying modern astronomical work, together with the results of these 
investigations. The facts are presented clearly and concisely, so that, in only 161 pages of 
textual material, the essential facts of astronomy are treated in a surprisingly comprehensive 
manner. The book has been prepared with care, so that the accuracy is high, even for a college 
text, and fourteen appendixes increase the value of the book for reference purposes. 

For a class in descriptive astronomy, composed primarily of nonscience majors, the reviewer 
would favor a fuller text, such as Duncan’s earlier text, Astronomy, even if certain sections of 
the book must be omitted. This type of student does better with a relatively full treatment of 
whatever points are studied. The fact that the arrangement of the new text requires treatment 
of a subject, such as precession, in two or more places, may increase the difficulty for these stu- 
dents. 

For a class composed primarily of majors in mathematics and science, especially a small 
class, the new text appears well suited. The arrangement is in agreement with scientific pro- 
cedure, helping the students to distinguish the facts of science, strictly speaking, on the one 
hand, from hypotheses and speculation, on the other. The brief treatment should cause these 
students no real difficulty, and the comprehensive coverage in a short space will be an advantage, 
especially in wartime. Professor Duncan has made a real contribution in writing a text for this 
type of course. 

C. C. WYLIE 


University of Towa 


National Geographic Society-National Bureau of Standards Solar Eclipse Expedition of 1940 to 
Brazil: Contributed Technical Papers of the National Geographic Society. (“Solar Eclipse 
Series,” No. 2.) Washington, 1942. Pp. 97. 

The Contributed Technical Papers of the National Geographic Society contain an account of 
the results secured in the joint expedition of the National Geographic Society and the National 
Bureau of Standards to observe the total eclipse of the sun on October 1, 1940. The station was 
Patos, Brazil. Although clouds partially obscured the total eclipse, much valuable information 
was obtained in spite of the meteorological handicaps. The publication consists of a foreword 
by Lyman J. Briggs, director of the National Bureau of Standards, followed by eight papers by 
various members of the expedition. 

In the first of these papers, Dr. Irvine C. Gardner discusses the design and construction of 
eclipse apparatus. The astronomical equipment consisted of two cameras for the corona, one 
fitted with polarizing apparatus for coronal studies, and two spectrographs, which together 
were ready to record the complete spectrum from 3000 to 10,000 A. The lenses of the corona 
cameras were coated with nonreflecting layers of calcium fluoride for the purpose of diminishing 
secondary reflections. 

The operation of the equipment was automatic, or at least semiautomatic, a feature that is, 
in the reviewer’s opinion, extremely desirable for eclipse equipment. Electric motors, fortu- 
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nately, have no tendency to get ‘‘buck fever’’ at the moment of totality. The corona cameras 
embody a number of excellent features. One wonders, however, whether the polaroid might 
not have been more advantageously placed just in front of the sensitive film rather than in 
front of the primary lens. 

In 1932 the reviewer designed a spectrograph in which the film was moved, semiautomati- 
cally, parallel to the dispersion between exposures. In 1936 he used a fully automatic camera, 
with wide film moved perpendicularly to the dispersion. The second method seemed definitely 
preferable, although the limitations produced by the finite width of the film and the conse- 
quently short spectrum were serious. Gardner has circumvented the difficulties very neatly 
with a camera in which the film moves diagonally, relative to the dispersion. As a result, he is 
able to record 40 inches of spectrum upon a film only 93 inches wide, with an advance of ap- 
proximately 12 inches between exposures. 

Dr. C. C. Kiess, of the National Bureau of Standards, in the next paper reports on the flash 
spectrum secured with the eclipse instruments. A fortunate break in the clouds near the end of 
totality enabled the scientists to obtain a record from about 3000 A to about 5000 A with a dis- 
persion of 2.4 A/mm in the first order. The spectra were taken a few seconds after the end of 
totality. Hence the flash lines appeared as crescent-shaped projections extending beyond the 
belt of continuous spectrum. Kiess has measured the wave lengths, estimated the flash intensi- 
ties, and recorded the character of the various lines. To designate the appearance, Kiess uses 
the letters a, 6, or c, according to the degree of line reversal in the strip of Fraunhofer con- 
tinuum. The classification is of undoubted astrophysical significance, but it seems unfortunate 
that Kiess has used a classification that resembles the one introduced by Menzel, without 
adopting the earlier classification. As a matter of fact, Kiess’s systenr is very similar to that of 
Menzel for types a and 6. But Kiess’s term c is generally denoted by d in the earlier classifica- 
tion. Kiess has failed entirely to note the distinctive behavior of the rare-earth lines which led 
Menzel to introduce an intermediate class c, which has no counterpart in Kiess’s system. Con- 
siderable confusion may result from the double classification. 

The next paper, by Stair and Coblentz, of the Bureau of Standards, describes the methods 
of calibration of the standard lens which was employed as a comparison for the flash spectrum. 
Kiess and Humphreys then discuss the determination of flash-spectrum intensities, microphoto- 
metrically, in what is unquestionably the most important paper of the series. They have in- 
vestigated only a few lines, H and K of Cau, the 6 lines of Mg 1, and the lines of the Balmer 
series of hydrogen from H8 to H¢. The authors find that for the metallic lines the well-known 
square-root law holds for the intensities. In view of the fact that their measures refer to very 
low-level flash, the result is not surprising. 

The authors have tabulated the hydrogen intensities in the customary manner, with 178 
fixed at unity, and have compared their measures with those of Cillié and Menzel on a similar 
scale. The intensities of the earlier members of the Balmer series, as given by Menzel and 
Cillié, are admittedly very uncertain, the values having been derived for the most part by 
extrapolation. 

The emission,! per frequency interval dv from a unit volume of gas in thermodynamic equi- 
librium, containing V; hydrogen ions and N, electrons per cm‘, is 


where gu is the Kramers-Gaunt factor, 7. the electron temperature, and n’ the quantum mem- 
ber of the lower level. Similarly, the total emission for a line involving the transition » to 7’ is 


(2rmk)3/2? mc8 3 V3" 


= 


Near the series limits g)~ gu, y™~ vn’ and n > ©. Hence the ratio of the intensity of a high- 
series line to the emission per unit frequency of the continuum at the limit is 


™ Menzel, Ap. J., 85, 330, 1937. 
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If, then, in any spectrum we arbitrarily set Exn’ = 1, for convenience, the high-series members 
that flow without break into the continuum will follow the law 


Baal = = 


n? 


The intensities of the remaining members should be referred to the continuum as a standard. 


TABLE 1 
| | OBSERVED INTENSITIES | THEORETICAL INYENSITIES 
Baker and Menzel 
| Live | Menzel and Cillié, 1932 
B A 
| K and H 
P, G, P, | = 10,000° | 10,000° 10,000° 
8 | 12.07* 11.90" 11.45 11.87} 12.03} 12.45 | 11.67 | 11.43 
y | 11.48* | 11.19 11.34) 11.75] 12.01 | 11.38 11.19 
6 | 11.14* | 11.08 10.84 10.98 11.52) 11.70 | 11.17 | 11 o1 
| 10.92* 1093 11.34) 11.47 | 11.00) 10.8 
| 10.73* | 10.89 | 11.14] 11.24 | 10:82 | 10/68 
7 | 10.66% 10.67 | 10.69] 10.50 | 10.70| 10.99] 11.07 | 10.69 | 10.56 
6 | 10-42% 10.45 | 10.45 | 10.46 | 10.66| 10.86] 10.92 | 10.57 | 10.45 
| | 10.39 | 10.38 | 10.38] 10.73} 10.78 | 10.45 | 10 33 
| 10.33* 10.36 | 10.34 10.29) 10.62} 10.66 | 10.37 | 10 23 
| 10.20* 10.31} 10.29 10.26} 10.52} 10:55 | 10.27| 10 14 
| 10.19% 10.21 | 10.24 10.17 10.43! 10.46 | 10.17! 10.07 
10.11 10.13 10.15) 10.33] 10.35 | 1007] 9.97 
: 10.03 | 10.07 9.99, 10.25! 16.27 | 10.00! 9 90 
17. 10.00 10.01; 9.98) 995) 10.18] 10.20 | 994 | 9.84 
r 991 9.95| 991 986) 10.09] 10.10 | 9 85 | 9.75 
p 989 9.89} 10.02} 10.03 | 9.79) 9 69 
| 982 9.79} 9.83) 9.85) 9.96] 9.07 | 9:73 | 9.64 
977| 9.79 979} 9 | 9.90 | 9.66} 9.57 
v | 971 | 9.70) 9.75} 9.72] 9.83} 9.84 | 9.61] 9.52 
| 971 | 965] 9.65} 970} 975! 976 | 9 54 | 9.45 
x | 9.62 9.50 ]....... 967! 967 | 938 
25... 9.49 9.59) 961) 9.61 | 9.39} 9.32 
w | 9.57 9.45|....... 9.55} 9.57] 9.57 | 9.35} 9.27 
| 9.51 9 40 9.52] 9.52 | 9.30] 90.22 
| 9 48 9.41, 947) 9:47 | 9:26 9:18 
30 9.39} 939 | 9 9.10 


| 
| 


The observational data for the 1932 eclipse clearly show that the chromosphere lines at 
n = 25-30 do grade properly into the continuum. Hence, setting R = 3.29 x 10, we may 
take 
log Ex, 9.39. 


This procedure is infinitely more satisfactory than the usual one of setting the intensityof HB 
equal to unity, although we can use it only if the measures include the continuum at the series 
limit. The reviewer presents in Table 1 a summary of the observational material with a com- 
parison of theoretical values adjusted to the above scheme. Columns 1 and 2 designate the 
Balmer lines by upper quantum number and customary notation. Columns 3, 4, and 5 sum- 
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marize the observations of Cillié and Menzel? from the 1932 eclipse. The designations P; and G, 
refer to two spectrograms taken simultaneously with prism and grating; P» refers to a prism 
film of a higher chromospheric level wherein the lines are weaker than in the other spectra. 
The intensities marked with an asterisk represent extrapolated values. Column 6 gives the 
values of Kiess and Humphreys, from the work under review. Column 7 lists intensities meas- 
ured by Thackeray* at the 1936 eclipse. 

Except for the first four or five lines of the series, the agreement is excellent. Over the range 
from H6@ to H¢ the mean discrepancy, taken without regard to sign, is only 0.03 in the logarithm. 
Also, since self-reversal may be present in varying amounts at the different eclipses and at 
different levels at the same eclipse and since its action would be greatest in the lower members, 
one has no right to question the relative accuracy of the different determinations. Changes of 
excitation conditions also have their greatest effect at Ha and HB. 

The remaining columns give values calculated from theory. The theoretical intensities for 
T = ©, from the Schroedinger-Sommerfeld formula, come next. Column 9 gives the intensities 
from an assembly in thermodynamic equilibrium at 10,000°. Columns 10 and 11 list the the- 
oretical intensities for two cases discussed by Baker and Menzel,‘ respectively, for complete ab- 
sorption of higher Lyman members (8B) and for no absorption of Lyman members, i.e., for a pure 
capture spectrum (A). 

The authors state that the Schroedinger-Sommerfeld formula for the line intensities of the 
hydrogen spectrum ‘‘applies to atoms undisturbed except for exposure to radiation of the Lyman 
series.”’ A more accurate statement would be that this formula applies to the emission of radia- 
tion by atoms whose population in the various levels is proportional to the weight of those 
levels; in other words, it applies strictly to the case of an infinite exeiting temperature. A cor- 
rection for the external condition must be introduced either by a Boltzmann factor or by the 
evaluation of the capture coefficients and mean lifetimes, as has been done by various theorists 
working on problems of planetary nebulae. The formula intensities, if uncorrected, have no 
physical significance in conjunction with the observations, except possibly near the series limit, 
where small differences in excitation potential have no marked efiect on the partition of atoms 
in the various levels over small ranges of energy and wave lengths. 

The reviewer is not in entire agreement with the astrophysical discussion of Kiess and 
Humphreys. It is true that the chromosphere presents some problems similar to those of the 
gaseous nebulae and emission-line stars. But the authors’ statement that “the hydrogen spec- 
trum in all these astronomical sources is regarded as essentially a recombination spectrum, 
little influenced by absorption of Lyman radiation from the central star or reabsorption of this 
radiation from the gas itself,’ is not confirmed by the observations. Also, it is by no means 
certain that the gas is too tenuous for collisions to play a significant part. The Balmer decre- 
ment, as measured by the ratio of the intensities of low- or intermediate-series members to the 
intensities of high-series members, is far too small to be explained on the basis of any Boltzmnan 
distribution, even at infinite temperature. As Menzel and Cillié pointed out long ago, there is an 
underpopulation of intermediate atomic levels which may amount to as much as a factor 
of 3 or 4. 

An analogous underpopulation is seen in the spectra of the planetary nebulae, where a very 
marked drop in the intensity of the continuum occurs where the lines overlap at the head of the 
series. For the chromosphere, the drop is by no means so pronounced, but the intensity meas- 
ures clearly prove its existence, i.e., the entries in the observed intensity columns tend to be 
less than those in columns 6 or 7. A theoretical investigation, however, shows that for a pure 
capture spectrum (col. 9) the drop is even more marked than observed for the planetary 
nebulae, whose intensities conform fairly well to the values tabulated in column 8. To account 
for the Balmer decrement, we shall have to suppose that reabsorption of the Lyman radiation 
is considerable both in the chromosphere and in the nebula. Since the lower members of even 
the Balmer series are probably partially reversed, the optical depth of the chromosphere in both 
the Lyman series and the continuum must be very high. The observed intensities are inter- 
mediate between those predicted on the capture theory, either with or without Lyman absorp- 
tion, and those predicted on the basis of pure thermodynamic equilibrium for infinite tempera- 
ture. Lower temperatures merely increase the decrement and tend to heighten the discrepancy. 


2 Harv. Cire., No. 410, 1935; Ap. J., 85, 88, 1937. 
3. M.N., 97, 672, 1937. ‘Ap. J., 88, 52, 1938. 
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It seems entirely reasonable to suppose that atomic collisions, despite the low densities, may 
play some role in redistributing atoms among the high atomic levels. Furthermore, self-absorp- 
tion of the lower members of the Balmer series, which is not taken into account in the elemen- 
tary theory, is in the right direction to account for any discrepancy. The fact that prominence 
and even chromospheric phenomena are clearly visible on the disk in Ha spectroheliograms is 
in itself the confirmation of the suggestion. Further theoretical work along these lines is clearly 
indicated. 

The advantage of the tabulation system here adopted is evident when we compare the the- 
oretical intensities with the observed. The Balmer decrement is not dissimilar for the two 
cases. But, if we follow our first impulse to add a constant to the figures tabulated in columns 
10 and 11, we must remember that we are also multiplying the predicted intensity of the con- 
tinuous spectrum by a factor. For the chromosphere this multiplication would predict a con- 
tinuum more intense than that observed. The reviewer reiterates the conclusion by Menzel and 
Cillié that the chromospheric spectrum is not predominantly one of pure capture and that con- 
siderable excitation by Lyman absorption or by electron impact must take place. 

In the next succeeding paper, Paul A. McNally, S.J., director of the Georgetown College 
Observatory, discusses the determination of contact times of the total solar eclipse from photo- 
graphs of the partial phases. In spite of the fact that the clouds interfered seriously, McNally 
found it possible to determine accurate contact times from the measured lengths of the oc- 
culted chords and the widths of the crescents. 

E. O. Hulburt, of the United States Naval Research Laboratory, in the following paper, 
discusses the determination of sky brightness at Patos, Brazil, during twilight and during the 
total solar eclipse. In the final scientific paper of the series, T. R. Gilliland, of the National Bu- 
reau of Standards, discusses the radio observations. One of the of the primary results derived 
from the observations was the determination of the order of magnitude of the recombination 
coefficients for the ions of the E, F;, and F, layers. Gilliland finds that there is some evidence of 
an ionizing force extending around the edge of the moon prior to reappearance of the sun after 
totality. The detection of this interesting phenomenon was made possible only because of the 
rapidity of the recording. The technique of observation well deserves close scrutiny by other 
prospective observers of eclipse-ionosphere phenomena. The last paper of the series is a narra- 
tive of the expedition, by Irvine C. Gardner. 

The complete set of papers represents a marked contribution. The generosity of the National 
Geographic Society in making the expedition possible has been amply rewarded by the excel- 
lence of the results. It is unfortunate indeed that the eclipse was partially cloudy; for the general 
value of the results obtained, despite the handicap, merely indicates how important the ob- 
servational data would have been had the full program been carried out. 


DonaLp H. MENZEL 
Harvard Observatory 


The Practical Essentials of Pre-training Navigation. By WILLIAM T. SKILLING and RoBErt S. 
RICHARDSON. New York: Henry Holt & Co., 1942. Pp. v+113. $0.75. 


The problem of supplying large numbers of flyers for the Army and Navy is a serious one in 
wartime, and the length of their training period is necessarily drastically curtailed, so that the 
student pilot gets his instruction under considerable pressure. The preflight-training courses 
which have been instituted in many high schools and colleges have the dual purpose of easing 
this strain to some extent and of providing a background to make the methods and procedures 
less sets of rules to be learned and applied. Begun in Britain in 1940, when the need for pilots 
was great, the plan has recently had considerable support in the United States. It is now well 
past the experimental stage and will continue to expand until the end of the war. Although no 
one can predict what the post-war world will be like, it is certain that aviation will play an im- 
portant part and that the high schools and colleges will be called upon to supply the physical 
and mathematical training for navigation, aerodynamics, and meteorology. 

It is the aim of the authors of this small volume to provide a text for certain of the subjects 
covered in such a course. The subtitle, Celestial Navigation—Meteorology—Map Projection, in- 
dicates the exact subjects treated. The reader is supposed to possess an almanac and a source of 
the correct time (not nearly so drastic an assumption as it would have been some years ago), 
and he may or may not have a sextant. In presenting the material from this standpoint, the 
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authors are able to restrict themselves to the general fundamental ideas involved and to avoid 
a great part of the complicated details of the use of instruments and equipment. The first 
chapters treat the spheres of the earth and the sky, the co-ordinates required to locate points 
upon them, and the different ways of attempting to represent a section of a sphere upon a plane 
surface. The chapter on ‘‘Time”’ is presented as clearly as such a puzzling subject can be, al- 
though it is not very likely that the practical navigator of the future will have much use for ap- 
parent time, sidereal time, the equation of time, or right ascension, and the marine navigator can 
even dispense with sidereal hour angle. 

The authors present a celestial fix as the intersection of two circles of position whose centers 
are the subsolar, lunar, or astral points and then show how to draw a line of position as a section 
of such a circle. It is to be hoped a treatment of this type should remove much of the mystery 
of celestial navigation. The last chapter is devoted to a short summary of the fundamentals of 
meteorology. 

Since the subject matter is specifically limited to celestial navigation, one cannot complain 
at the absence of any mention of plane or great-circle sailing, wind drift, or the plotting of ter- 
restrial lines of position, which would have to be included in a complete pretraining course in 
navigation. 

The illustrations are very carefully and clearly drawn and should contribute to the reader’s 
understanding, although he will be puzzled by the authors’ method of constructing a mercator 
projection on page 12. 

The instructor of a high-school course in navigation will find in this volume a manual for his 
students which presents the fundamentals of the subjects clearly and avoids the confusing de- 
tails which make the more complete texts difficult for the beginner... 

Joun Titus 
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